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FORWARD 

The  work  reported  herein  summarizes  the  final  research  activity 
in  the  Beach  Environment  Task  of  a  program  at  ERIM  entitled  "Coastal 
Remote  Sensing  Investigations."  The  work  has  been  performed  under  the 
guidance  of  Mr.  Hans  Dolezalek  of  the  Office  of  Naval  Research,  Code 
460. 


This  study  was  carried  out  in  the  Applications  Division  of  the 
Environmental  Research  Institute  of  Michigan  (ERIM).  The  principal 
investigator  for  the  program  was  Dr.  David  R.  Lyzenga  while  Mr.  Eric 
S.  Kasischke  served  as  task  leader  for  the  project.  Mr.  Robert  A. 
Shuchman  supplied  technical  expertise  for  running  the  MOGS  algorithm. 
Dr.  Gwynn  Suits  contribution  Section  6  dealing  with  the  extension  of 
the  MOGS  algorithm  to  silts  and  clays. 
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INTRODUCTION 


This  report  summarizes  the  final  research  activity  in  the  Beach 
Environment  task  of  a  multi-year  study  at  ERIM  concerned  with  the 
problem  of  Coastal  Remote  Sensing.  The  purpose  of  the  Beach 
Environment  task  was  to  study  ways  in  which  beach  environmental 
parameters  (soil  particle  size,  moisture  content,  and  mineralogical 
composition)  could  be  remotely  sensed.  The  remotely  sensed  beach 
parameters  are  then  used  as  inputs  in  a  trafflcabil ity/mobil ity  model. 

Included  in  this  report  is  a  summary  of  previous  work  accom¬ 
plished  under  this  program,  as  well  as  specifically  reporting  on  final 
year  activities.  Scientific  findings  from  this  ONR  supported  study 
have  been  published  in  the  open  literature.  The  appendix  to  this 
report  Includes  these  reprints  as  background  information. 

In  the  final  year  of  this  effort,  several  analyses  were 
performed.  A  more  accurate  and  environmentally  applicable  mineralogy, 
moisture  and  grain  size  algorithm  (MOGS)  was  developed  using  actual 
sand  samples  and  inputs  from  the  AQUASAND  model.  The  MOGS  algorithm 
first  classifies  the  sand  sample  into  one  of  five  basic  beach 
classifications.  These  classifications,  based  on  mineralogical 
breakdowns,  were  performed  via  a  minimum  vector  length  technique. 
Next,  the  MOGS  algorithm  estimates  the  moisture  and  grain  size  using 
equations  based  on  actual  and  AQUASAND  data. 

The  MOGS  algorithm  was  evaluated  using  "unknown"  beach  sand 
samples.  The  MOGS  algorithm  correctly  identified  the  beach  mineralogy 
in  3  out  of  6  cases.  The  percent  moisture  estimates  predicted  by  the 
algorithm  were  in  good  agreement  with  the  actual  moistures.  There 
seems  to  be  some  doubt  as  to  whether  the  MOGS  algorithms  can  be  used 
to  accurately  predict  grain  size.  This  failure  is  not  serious  from  a 
beach  trafficabil ity  standpoint,  however,  because  recent  studies  have 
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shown  that  beach  trafficabil Ity  models  are  rathar  insensitive  to  grain 
size,  when  all  other  factors  are  considered  (Thomson,  et  al.,  1980). 

The  MOGS  algorithm  was  then  modified  so  that  it  could  be  used  to 
classify  data  collected  b>  the  ERIM  airborne  multi  spectral  scanner. 
Data  from  two  sites  were  processed.  One  site  was  on  the  shores  of 
Lake  Michigan,  near  Muskegon,  Michigan,  while  the  other  site  was  at 
Panama  City,  Florida. 

The  final  sections  of  this  report  discuss  the  use  of  satellite 
data  as  inputs  into  the  algorithms  and  the  development  of  a  MOGS  type 
algorithm  for  soils  other  than  sand. 
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BACKGROUND 

The  ability  to  remotely  estimate  composition  of  surface  terres¬ 
trial  materials  has  been  demonstrated  by  many  workers,  and  these 
results  may  have  direct  transfer  value  to  beaches,  where  the  composi¬ 
tion  is  similar  to  that  studied  by  geologists.  The  most  prominent 
reports  in  this  area  have  been  by  Vincent  (1972,  1973),  Vincent  and 
Thomson  (1972),  and  Vincent  et  al .  (1972).  These  authors  have 
explained  the  physical  and  optical  properties  of  various  rock  and  soil 
types  as  well  as  demonstrated  the  capabilities  of  a  thermal  infrared 
scanner  in  accurately  classifying  silica-containing  soil  and  rock  out- 
croppi/igs  and  areas  with  surface  iron  stain.  Discrimination  is  gener¬ 
ally  based  on  either  atomic  or  ionic  absorption  in  the  shorter  wave¬ 
lengths  (0.35  to  2.5  um)  and  by  lattice  absorption  at  longer  wave¬ 
lengths  (8  to  14  um  reststrahlen  bands  in  silicate  and  carbonate 
minerals).  These  effect  the  spectral  reflectance  or  emittance  of 
materials. 

Reststrahlen  spectral  emissivity  variations  with  surface  moisture 
have  been  studied  by  Watson  (1970).  He  determined  that  as  a  thickness 
of  water  increased  from  0  to  30  um  on  a  polished  quartz  slab,  emit¬ 
tance  increased  at  the  principal  quartz  reststrahlen  wavelengths  of 
8.5,  9.0,  and  12.5  um.  This  may  be  an  important  way  to  detect  sample 
moisture  content,  and  the  effect  will  almost  certainly  influence  com¬ 
positional  maps  made  of  wet  areas  using  the  techniques  reported  by 
Vincent  (1972,  1973). 

Another  qualitative  observation  is  that  the  size  of  silicate 
mineral  grains  affects  the  reststrahlen  emissivity  features.  Hovis 
and  Callahan  (1966)  have  demonstrated  qualitatively  that  when  the 
grain  size  of  a  sample  is  reduced  from  a  solid  sample  to  particles 
with  diameters  of  1-2  mm,  0.105-0.250  mm,  and  less  than  0.038  mm,  the 
measured  reflectance  progressively  increases.  Hunt  and  Vincent  (1968) 
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have  attempted  to  explain  this  phenomena  in  terms  of  a  specular 
surface  reflectance  component  (R$)  and  a  reflectance  component  which 
is  due  to  radiation  that  has  been  transmitted  through  part  or  parts  of 
the  sample  before  re-emerging  (Ry).  The  total  sample  reflectance  is 
the  sum  of  R$  and  Ry.  Further,  the  ratio  Rg/Rv  depends  on  grain  size, 
because  R$  and  Ry  depend  differently  on  the  absorption  coefficient,  a. 
Thus,  the  reflectance  or  emittance  of  the  sample  will  depend  on  grain 
size. 


Because  of  the  promising  research  described  above,  ERIM  began,  in 
1974,  to  relate  .grain  size,  degree  of  sorting,  composition,  and  water 
content  of  beach  sands  to  the  spectral  properties  in  the  visible, 
reflective  infrared,  and  thermal  infrared  regions.  The  objective  was 
to  find  quantitative  relationships  between  the  spectral  reflectances 
of  beaches  and  their  physical  properties.  During  the  first  two  years 
of  the  program,  50  beach  sand  samples  from  a  variety  of  environments 
were  collected.  These  samples  were  analyzed  in  terms  of  mineralogy, 
grain  size,  degree  of  sorting  and  degree  of  iron  staining.  In 
addition  the  spectral  reflectance  of  each  was  determined  in  the  0.35 
to  2.5  um  range  using  the  Cary  14  spectrophotometer.  A  number  of 
different  moisture  contents  were  used  to  simulate  different  regions  in 
the  land/water  interface.  Thermal  infrared  measurements  were  made  in 
the  8-14  urn  band  on  some  of  the  sand  types,  however,  further  analysis 
was  restricted  to  visible  and  near  IR  wavelength  because  of  moisture 
attenuation  in  this  region. 

Several  regression  equations  were  developed  at  the  end  of  the 

second  year  in  an  effort  to  predict  grain  size  and  moisture  from 

spectral  data.  Among  them,  a  two  channel  moisture  algorithm  was 

2 

developed  which  resulted  in  a  coefficient  of  variation,  R  of  0.79  and 
a  standard  error  of  4.5%.  Another  moisture  algorithm  that  used  five 
spectral  channel  ratios  was  capable  of  predicting  moisture  content 
with  a  standard  error  of  4.0%  and  a  Rc  of  0.86.  For  grain  size,  one 
algorithm  using  10  ratios  was  found  capable  of  predicting 
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the  Wentworth  grade  scale  of  beach  sands  with  an  R  of  0.83.  Although 
the  results  appeared  very  good,  the  physical  reasoning  for  the 
spectral  channel  selection  and  success  of  these  algorithms  was  felt  to 
be  somewhat  in  doubt.  In  addition  the  large  number  of  factors  in  the 
regression  coupled  with  the  relatively  small  data  set  reduced  the 
significance  of  the  results.  These  equations  did,  however,  show  the 
feasibility  of  developing  an  algorithm  which  could  successfully 
predict  moisture  and  grain  size  from  spectral  data. 

In  order  to  better  control  the  variable  parameters  of  beach  sands 
and  to  more  economically  create  the  necessary  data  sets,  work  began  in 
the  third  year  on  an  adaptation  of  the  Suits  radiative  transfer 

vegetative  canopy  model  to  the  beach  sand  situation.  The  initial 
model  ,  SANDREF,  was  verified  by  the  end  of  the  third  year  on  two  dry 
sands;  a  carbonate  type  and  a  pure  quartz  type.  In  the  process  of 
developing  this  unique  modeling  approach,  the  scattering  coefficient 
(s),  forward  scattering  (FS),  and  absorption  coefficient  (a)  of 
sixteen  common  beach  forming  minerals  were  derived.  This  was  the 
first  time  such  a  measurement  had  been  achieved  for  these  minerals. 
In  order  to  develop  these  parameters  three  thin  sections  of  each 
mineral  of  interest  had  to  be  cut.  Each  section  had  to  be  slightly 
thinner  than  the  one  before  it.  Using  the  reflectance  and 
transmittance  of  each  of  the  three  thin  sections,  an  iterative  curve 
fitting  procedure  was  used  to  develop  the  a,  FS,  and  s  parameters. 

These  parameters  were  then  used  to  predict  the  reflectance  and 

transmittance  of  a  given  mineral  at  any  desired  thickness.  The 

reflectance  and  transmittances  were  then  entered  into  the  SANDREF 
model . 

After  SANDREF  was  proven  to  work,  it  was  modified  to  utilize 
moisture  input  to  simulate  the  spectra  of  sand  which  contains  water. 
The  new  program,  called  AQUASANDREF,  included  a  moisture  algorithm. 
The  two  modeling  programs  were  then  combined  at  the  beginning  of  the 
fourth  year  into  a  single  all-purpose  program  called  AQUASAND 
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(Shuchman,  et  al.,  1978).  This  composite  program  can  theoretically 
handle  any  beach  sand  situation.  It  should  be  stated  that  the 
AQUASANO  model  is  not  designed  to  predict  sand  parameters  from 
reflectance  spectra  but  rather  to  generate  spectra  from  known 
parameters.  These  spectra,  In  conjunction  with  measured  spectra  may 
be  used  to  generate  a  predictive  algorithm. 

Specifically,  during  year  four  a  new  method  of  modeling  the  iron 
stain  frequently  found  as  a  coating  on  sand  grains  was  developed,  the 
two  models,  AQUASANDREF  and  SANDREF,  were  combined  and  shortened  into 
a  much  more  efficient  form  (AQUASAND)  for  beach  sand  modeling.  This 
version  handled  modeling  for  both  wet  and  dry  sands,  and  a  total  of 
seven  diverse  beach  types  were  modeled  successfully  using  separate 
mineralogical  and  moisture  inputs  to  the  AQUASAND  model.  This  work 
leads  into  the  research  reported  in  this  report:  the  development  of  a 
mineralogy,  moisture  and  grain  size  algorithm,  and  the  verification  of 
this  model  on  laboratory  and  aircraft  produced  data. 
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MOGS  ALGORITHM 

Through  the  use  of  the  AQUASAND  model,  we  were  able  to  generate 
the  spectral  signature  of  beach  sands  by  varying  the  mineralogy, 
moisture  and  grain  size.  The  goal  of  this  study  was  to  reverse  the 
AQUASAND  model,  i.e.,  given  the  reflectance  spectrum  of  a  sand  sample, 
generate  the  mineralogy,  grain  size,  and  moisture  of  that  sample. 
This  section  will  discuss  the  development  of  such  an  algorithm. 

Examples  of  actual  versus  AQUASAND  generated  spectra  are  given  in 
Figure  1.  From  inspection  of  these  examples,  one  can  see  the  ability 
of  the  AQUASAND  model  to  closely  predict  the  spectra  of  certain  beach 
sands.  Further  examination  of  the  beach  sand  spectra  revealed  that 
the  mineralogy  had  by  far  the  greatest  influence  on  bulk  reflectance. 
So  great  is  this  influence  that  it  tends  to  mask  the  more  subtle 
features  of  changes  in  grain  size  and,  to  a  much  lesser  extent, 
moisture.  It  was  decided  that  in  order  for  the  algorithm  to  handle  a 
broad  variety  of  sand  types  and  still  maintain  the  resolution  needed 
to  detect  small  spectral  features,  a  preprocessing  classification  of 
mineralogy  would  be  necessary. 

In  order  to  discriminate  mineralogy,  a  vector  length  decision 
framework  was  used.  The  concept  is  developed  as  follows. 

Suppose  that  there  are  two  points,  A  and  B,  located  in  two 
dimensional  space.  The  distance,  or  vector  length,  L,  from  A  and  B 
can  be  expressed  in  terms  of  the  X  and  Y  coordinate  locations  of 
points  A  and  B  as: 


L  ■  (XA  -  XB)2  ♦  (Ya  -  Y,)*. 


0) 
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FIGURE  1 •  COMPARISON  OF  MODEL  GENERATED 
AND  SPECTRALLY  MEASURED  REFLECTANCES  FOR 
SEVEN  COMMON  BEACH  TYPES 
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This  is,  of  course,  related  to  the  Pythagorean  Theorem.  Now  suppose 
we  have  a  p-dimensional  system  with  A  and  B  located  in  each  dimension, 
the  vector  length  can  be  expressed  as 


where  is  the  location  of  point  A  in  the  ith  dimension  and  is 
the  location  of  point  B  in  the  ith  dimension. 

This  rationale  can  be  used  to.  classify  some  point,  T,  as  being 
the  member  of  one  of  n  classes  (A^,  j  =  l,n),  by  finding  the  minimum 
vector  length  from  T  to  A.(j  =  l,n).  In  other  words  T  is  said  to  be  a 

J 

member  of  the  class  which  is  closest  to  it,  on  the  average  across  all 
p  dimensions.  The  minimum  vector  length  is  defined  as 


Notice  that  equation  3  has  no  provision  for  variability  in  the  n 
classes,  therefore,  Lm^n  is  chosen  as  being  the  shortest  linear  vector 
length.  If  each  class  has  the  same  variability  associated  with  it, 
this  causes  no  difficulty.  In  this  experiment,  however,  there  were 
considerable  differences  in  variability  between  the  classes  so  that  a 
modification  of  equation  3  had  to  be  made.  The  standard  deviation 
(SD)  was  used  to  modify  the  distance  between  T  and  Aj  related  to  each 
dimension  thus  removing  the  effects  of  variability  from  each  class. 
This  normalized  minimum  distance  equation  is  expressed  as 


(4) 
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where  SD^  is  the  standard  deviation  associated  with  the  jth  class  in 
the  ith  dimension. 

In  the  application  of  this  method  to  the  classification  of 

mineralogy,  the  "dimensions"  are  spectral  bands  or  ratios  of  spectral 
bands  and  the  "classes"  are  mineralogical  types.  Eight  spectral  bands 
(Table  1)  and  all  possible  unique  ratios  of  those  spectral  bands  were 
used  to  classify  the  mineralogical  type  of  an  input  sand  as  cne  of 
five  categories  (Table  2).  The  object  was  to  make  each  category  as 
homogeneous  as  possible  so  that  the  moisture  and  grain  size 
regressions  which  followed  would  be  sensitive  to  small  scale  spectral 
changes. 

Using  the  AQUASAND  model  we  found  that  information  related  to 

moisture  content  of  sands  is  best  derived  from  the  spectral  region 
beyond  1.0  pm.  This  is  due  to  the  fact  that  the  spectral  reflectance 
of  sand  in  this  region  is  reduced  by  absorption  in  proportion  to  the 
amount  of  water  present.  Exceptionally  high  spectral  absorption  is 
noted  near  1.4  and  1.9  urn.  Although  the  spectral  reflectance  in  these 
regions  is  highly  correlated  to  moisture  we  did  not  consider  them 
since  atmospheric  absorption  prohibits  their  use  by  an  airborne 
sensor. 

Changes  in  grain  size  seem  to  manifest  themselves  most  clearly  in 
the  shorter  wavelengths  (.4-. 7  uni).  Grain  size  information  is  gained 
by  light  being  reflected  from  sand  grains  below  the  surface  through 
surface  grains.  The  transmittance  through  the  surface  grains  is 

reduced  by  internal  scattering  and  absorption  of  the  particle.  Since 
both  of  these  factors  are  dependent  on  thickness,  the  bulk  reflectance 
of  a  sand  is  dependent  to  some  degree  o-i  the  grain  size. 
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Table 


Table 


1.  The  8  spectral  bands  used  in  the  breakdown  of 
beach  mineralogy  into  1  of  5  categories.  In 
addition  to  these  8  bands  all  unique  ratio 
combinations  were  also  used. 


Band  # 

1 

2 

3 

4 

5 

6 

7 

8 


Wavelength  Range  (im) 

.43-. 47 
.47-. 49 
.51-. 53 
.53-. 56 
.59-. 63 
.80-. 90 
.90-1 .0 
1  .0-1  .1 


2.  The  5  potential  mineralogical  classifications 


Class  # 

1 

2 

3 

4 

5 


Description 

Iron  stained  Atlantic  coast  type 
Iron  stained  Michigan  coast  type 
Iron  stained  pure  quartz  type 
Heavy  mineral  type 
Carbonate  type 
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Theoretically  a  large  grain  sand  should  have  a  lower  reflectance  than 
a  small  grain  sand  of  the  same  mineralogical  composition  and  with 
similar  moisture  content.  According  to  our  measurements  this  appears 
to  be  the  case. 

This  grain  size  phenomenon  can  be  confounded  in  two  ways.  First 
if  there  is  no  scattering  or  absorption  within  the  grains  (i.e.,  a 
perfectly  clear  material  at  all  wavelengths)  there  can  be  no 
attenuation.  Fortunately,  even  in  our  purest  quartz  sands  there  were 
enough  impurities  and  inclusions  to  give  some  attenuation.  Second, 
the  sand  grains  may  be  opaque  and  thus  attenuate  too  much  light.  This 
appears  to  be  the  case  in  the  heavy  mineral  and  carbonate  beaches. 
Most  of  the  bulk  reflectance  for  these  two  types  was  due  to  surface 
reflectance  and  little  if  any  was  due  to  light  transmitted  through  the 
surface  grains  from  below.  We  were  unable  to  create  accurate  grain 
size  equations  for  these  types. 

Utilizing  the  physical  phenomena  discussed  above  we  were  able  to 
develop  multiple  linear  regression  equations  for  predicting  moisture 
in  all  five  mineralogical  classes  and  grain  size  for  three  of  the  five 
mineralogical  classes.  The  basis  for  all  the  regressions  except  one 
was  the  sample  group  corresponding  to  a  given  mineralogical  class. 
The  single  exception  was  the  grain  size  equation  corresponding  to  a 
pure  quartz  beach.  Our  samples  within  this  type  consisted  of  a  single 
grain  size  (0.22  mm)  and,  as  such  did  not  provide  an  adequate  basis 
for  regression  equations.  For  this  case  we  used  AQUASAND  generated 
spectra  to  simulate  a  wide  range  of  grain  sizes  in  order  to  add  grain 
size  variability  to  the  data  set. 

Seventeen  spectral  bands  between  0.4  and  2.5  urn  were  chosen  for 
use  in  the  regressions  (Table  3).  Within  the  17  bands,  only  those 
which  were  predicted  by  the  AQUASAND  model  to  be  most  informative  were 
used.  In  this  way  we  could  be  reasonably  certain  that  the  regression 
equations  would  respond  the  correct  parameter  and  thus  yield  accurate 
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Table  3.  The  17  bands  used  in  the  development  of  moisturs 
and  grain  size  regression  equations 


Band  # 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


Wavelength  Range  (urn) 

0.43-0.47 
0.47-0.49 
0.49-0.51 
0.51-0.53 
0.53-0.56 
0.56-0.59 
0.59-0.63 
0.63-0.67 
0.70-0.75 
0.75-0.80 
0.80-0.90 
0.90-1.00 
1 .00-1 .10 
1 .10-1 .20 
1 .20-1 .35 
1 .50-1 .85 
2.10-2.50 
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predictions.  The  predictive  equations  together  with  the  associated 
standard  errors  (SE),  and  coefficient  of  variation  (R  )  are  given  in 
Table  4. 


In  summary  the  MOGS  algorithm  (Figure  2)  represents  a  computer 
controlled  package  of  equations.  The  input  is  a  set  of  17  spectral 
reflectance  bands  obtained  from  an  unknown  sand.  Based  on  these 
bands,  the  sand  is  classified  as  being  a  member  one  of  five 
mineralogica'i  types.  Depending  on  the  mineralogical  type,  the 
appropriate  moisture  and  grain  size  (where  applicable)  equations  are 
applied  to  the  data.  The  output  from  the  MOGS  algorithms  is  the 
predicted  mineralogical  class,  the  predicted  moisture,  and  the 
predicted  grain  size. 
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Table  4.  Multiple  linear  regression  equations  for  the  prediction 
of  moisture  and  grain  size.  The  equations  are  listed  by 
mineralogical  class.  Grain  size  is  in  mm. 

Iron  stained  quartz  -  Atlantic  coast 

Predicted  moisture  %  *  67.964  -  65.046  j^.) 

S.E.  =  3 . 08% ,  R2  =  0.888 

Predicted  grain  size  =  6.87  -  3.4634  (Band  7)^  +  .0300  (Band  1)  + 

.01672  (Band  15) 

S  E.  =  0.13  mm,  R2  *  0.603 

Iron  stained  quartz  -  Michigan  coast 

Predicted  moisture  %  =  60.149  -  49.961  (ff^-jj)  -  2.226 

S.E.  =  2.56% ,  R2  =  .970 

Predicted  grain  size  *  0.6405  -  0.0152  (Band  5)  -  .0047  (Band  17) 

S.E.  =■  0.055  m,  R2  =  0.558 
Non-Iron  stained  quartz 

Predicted  moisture  %  *  127.02  -  65.159  -  64.054 

S.E.  =  2.12%,  R2  =  0.971 

Predicted  grain  size  =  1.158  -  2.328  (Band  10)  +  .3201  (j^jj-y)  + 

0.2858  (Band  16) 

2 

S.E.  and  R  not  applicable 
Carbonate 

Predicted  moisture  %  =  596.28  -  642  -  1.081  (Band  14)  + 

0.1538  (Band  17) 

S.E.  =  4.09%,  R2  =  .879  No  grain  size  equation 
Heavy  mineral 

Predicted  moisture  %  =  19.284  +  11.194  -  1.081  (Band  14)  + 

0.1538  (Band  17) 

? 

S.E.  =  4.09%,  R  =  .879  No  grain  size  equation 
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FIGURE  2.  FLOW  OIAGRAM  OF  THE  MOGS  ALGORITHM. 
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MOGS  ALGORITHM  TESTS 

It  is  recognized  that  the  MOGS  algorithm  was  generated  from  a 
very  limited  data  set.  It  is,  therefore,  necessary  to  test  this  model 
before  a  decision  is  made  on  its  usefulness.  Not  surprisingly,  when 
the  sand  samples  from  which  the  algorithm  was  developed  were  used  to 
test  the  model,  very  good  results  were  obtained.  This,  however,  was 
expected.  To  realistically  evaluate  the  model,  more  rigorous  tests 
had  to  be  conducted. 

Not  only  should  other  laboratory  generated  sand  spectra  be 
tested,  but  the  use  of  this  algorithm  on  multi  spectral  scanner 
collected  data  also  needed  to  be  evaluated.  This  section  deals  with 
two  such  tests. 

4.1  LABORATORY  TESTS 

Six  unknown  (to  the  technician  conducting  the  tests)  sands  were 
collected  for  use  in  testing  the  MOGS  algorithm  under  laboratory 
conditions.  Three  of  these  sands  were  from  Atlantic  Coast  quartz 
beaches,  two  were  from  carbonate  beaches  and  one  was  from  a  Michigan 
quartz  beach.  Each  sample  was  subdivided  into  subsamples.  These 
subsamples  were  mixed  with  varying,  known  amounts  of  water  to  create 
different  moisture  conditions. 

The  spectra  of  these  subsamples  were  next  measured  using  a 
Beckman  spectrophotometer.  The  output  from  the  Beckman  were  in  analog 
form;  therefore,  the  data  had  to  be  digitized.  Considerable  software 
development  was  necessary  to  transform  the  digitized  data  into  a  form 
which  was  compatible  for  input  into  the  MOGS  algorithm. 

Table  5  summarizes  the  results  for  the  sand  samples  analyzed  in 
the  laboratory.  Given  in  this  table  are  the  sample  number,  the  sand 
identification  (actual  and  predicted)  the  percent  moisture  (actual  and 
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predicted)  and  the  grain  size  {actual  and  predicted).  We  can  see  that 
the  MOGS  algorithm  correctly  identified  the  mineralogy  only  50  percent 
of  the  time.  The  reason  for  this  is  two-fold.  First,  the  MOGS 
algorithm  was  developed  using  a  very  limited  number  of  beach  sand 
samples.  Therefore,  when  confronted  with  a  sand  which  has  a  very 
different  composition  from  those  used  in  the  development  of  the  model, 
the  MOGS  algorithm  has  difficulty  classifying  that  sand.  Secondly, 
the  instrument  which  originally  measured  the  sand  samples  for 
development  of  the  MOGS  algorithm  was  a  Cary-14  spectrophotometer. 
The  instrument  used  to  measure  the  test  samples  was  a  Beckman 
spectrophotometer.  While  these  two  instruments  operate  using  the  same 
principle,  and  produce  quite  similar  results,  it  is  possible  that 
these  instruments  would  give  two  different  reflectances  for  the  same 
sample  over  very  narrow  wavelengths.  We  can  see  from  Table  1  that 
many  of  the  bands  used  in  the  mineralogy  classification  are  narrow 
bands;  a  small  variation  in  the  location  of  these  bands  between  the 
two  instruments  could  result  in  a  misclassification. 

The  moisture  results  were  quite  good.  The  correlation  of 
predicted  to  actual  moisture  is  0.93  (significant  at  the  0.01  level). 
The  moisture  predictions,  however,  appear  suspect  at  higher  moistures. 

The  grain  size  predictions  were  quite  poor  for  the  laboratory 
examples.  This  is  in  part  due  to  the  MOGS  algorithm's  misclassifica¬ 
tion  of  the  mineralogy  of  half  the  examples,  and  in  part  to  the 
limited  number  of  samples  used  in  developing  the  algorithm. 

4.2  MULTI SPECTRAL  SCANNER  TESTS 

The  MOGS  algorithm  was  tested  on  two  separate  sets  of  ERIM 
collected  multi  spectral  scanner  data.  One  set,  collected  at  Panama 
City,  Florida,  was  obtained  before  the  development  of  the  MOGS 
algorithm,  and  only  contained  enough  spectral  bands  of  data  to  provide 
a  moisture  map  of  the  test  area.  Data  collected  over  a  Lake  Michigan 
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test  site  was  obtained  during  the  development  of  the  MOGS  algorithm 
and  contained  the  proper  spectral  bands  for  the  MOGS  algorithm. 

4.2.1  PANAMA  CITY  TEST  SITE 

Multi  spectral  scanner  (MSS)  data  from  the  morning  of  23  April 
1976  was  interpreted  to  be  the  best  available  from  the  Panama  City 
test  site  for  input  into  the  MOGS  algorithm.  Unfortunately,  none  of 
the  Panama  City  collected  data  contained  enough  spectral  bands  for 
grain  size  analysis.  Additionally,  of  the  three  spectral  bands  needed 
for  the  moisture  equation  (see  Table  4)  only  two  were  available,  and 
thus  the  moisture  algorithm  had  to  be  modified  accordingly. 

Figure  3  shows  the  moisture  map  for  the  Biltmore  Beach  test  site. 
The  Biltmore  Beach  is  a  typical  coastal  beach  in  that  it  is  well 
drained,  and  on  a  sunny  day  like  the  one  when  the  data  was  collected, 
the  surface  layer  of  sand  is  virtually  moisture  free. 

4.2.2  LAKE  MICHIGAN  TEST  SITE 

The  ERIM  MSS  was  flown  over  a  portion  of  the  Lake  Michigan  shore¬ 
line  near  Pentwater,  Michigan  on  1  November  1978.  ERIM  scientists 
were  present  during  this  time,  collected  samples  were  for  laboratory 
analysis  of  grain  size  and  moisture  content.  The  MOGS  algorithms  was 
modified  to  conform  to  the  12  band  configuration  of  the  ERIM  scanner 
used  during  this  data  acquisition. 

The  actual  and  MOGS  predicted  moisture  and  grain  size  data  for 
the  field  collected  samples  are  presented  in  Table  6.  The  agreement 
between  the  actual  and  the  MOGS  predicted  parameters  appear  to  be 
quite  good,  but  the  moisture  prediction  appears  particularly  poor  for 
high  moisture  contents.  This  may  be  due  to  the  fact  that  wet  sand  at 
the  test  sites  appeared  to  exhibit  bi-directional  dependencies  (i.e., 
a  failure  to  behave  in  a  Lambertian  manner).  These  bi-directional 
characteristics  are  enhanced  at  low  sun  angles  and  are  not  accounted 


Table  6.  Comparison  of  actual  parameters  to  predicted 
classifications  by  the  MOGS  algorithm.  The 
sand  spectra  used  here  were  collected  by  the 
ERIM  MSS 


Sample 

Mineralogy 

Moisture  % 

Grain 

Size  (mm) 

Actual 

Predicted 

Actual 

Predicted 

Actual 

Predicted 

MICH1 

MICH 

22.3 

12.1 

.25 

.32 

MICH2 

MICH 

1 .0 

0.0 

.23 

.28 

MICH  3 

MICH 

8.0 

10.1 

.22 

.25 

MICH4 

MICH 

28.0 

16.9 

.26 

.35 

MICH5 

MICH 

5.0 

2.2 

.24 

.26 
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for  by  the  MOGS  algorithm.  Although  the  flight  took  place  at  1:30  EST 
the  sun  was  only  38°  above  the  horizon  on  November  1.  To  minimize 
bi-directional  reflectance  variations,  future  aircraft  flights  should 
be  made  during  complete  mid-altitude  (3000  m)  cloud  cover  or  sunny 
skies  with  the  sun  close  to  the  zenith  (summer  sun). 

Applying  the  same  moisture  and  grain  size  equations  used  in  the 
previous  analysis,  the  entire  Pentwafer  State  Park  beach  on  a  pixel  by 
pixel  basis  (in  this  case  1.5  x  1.5  meters)  was  classified  with 
respect  to  grain  size  and  moisture  content.  The  two  MOGS  generated 
digital  maps  (see  Figure  4)  show  the  predicted  moisture  and  grain  size 
distributions  on  the  beach  at  Pentwater.  The  ground  truth 
measurements  taken  at  the  time  of  flight  correlate  well  with  these 
images.  Figure  4  helps  to  demonstrate  how  an  entire  sandy  coastline 
could  be  analyzed  in  respect  to  moisture,  grain  size,  and  gross 
mineralogy  using  a  small  subsection  as  calibration. 
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MOISTURE  DISTRIBUTION  IMAGE  -  I. IGHTF.K  SHADE  INDICATES  HIGHER  MOISTURE 

CONTENT 


FIGURE  4 


0-57. 


5-207. 


20+7. 


GRAIN  SIZE  DISTRIBUTION  LIACi.  -  LIGHTER  SHADE  INDICATES  LARGER  GRAIN  SIZE 


.15-. 25  mm 


.  25-  .  40  mm 


40  mm 


PANCHROMATIC  AERIAL  PHOTOGRAPH 


DIGITAL  IMAGERY,  GENERATED  BY  THE  MOGS  ALGORITHM  SHOWING  THE  DISTRIBUTION 
OF  MOISTURE  AND  GRAIN  SIZE  ON  PENTWATER  BEACH  (PENTWATER  STATE  PARK, 
MICHIGAN).  WHITE  AREAS  ARE  EITHER  UNCLASSIFIED  REGIONS  OR  OPEN  WATER. 
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USE  OF  SATELLITE  MULTI  SPECTRAL  SCANNER  DATA 
AS  MOGS  ALGORITHM  INPUTS 

Based  on  the  parameters  of  current  spaceborne  (MSS)  systems  and 
the  specific  requirements  of  the  MOGS  algorithm,  there  appears  to  be 
very  limited  potential  in  using  satellite  data  as  inputs  into  the 
algorithm.  This  is  based  on  two  observations:  (1)  the  present  number 
of  spectral  bands  available  on  satellite  multi  spectral  scanners  are 
limited  to  three  or  four,  where  the  MOGS  algorithm  currently  needs  a 
large  number  of  spectral  bands  to  operate;  and  (2)  the  resolution  of 
spaceborne  scanner  systems  are  probably  too  coarse  to  image  beach 
environments. 

The  MOGS  algorithm  was  derived  using  eight  bands  of  information 
for  the  initial  mineralogical  classification.  After  this,  up  to  six 
more  bands  were  needed  to  produce  the  moisture  and  grain  size 
predictions.  Current  satellite  systems  (i.e.,  the  Landsat  series) 
have  only  four  channels  of  information  (see  Table  7).  From  Table  7, 
it  can  be  seen  that  future  proposed  satellite  systems  will  have  more 
channels  than  the  existing  Landsat  series.  But  even  the  six  channels 
of  data  which  both  Landsat-D  or  CZCS  will  have  do  not  come  close  to 
matching  the  14  channels  of  data  the  MOGS  algorithm  needs.  Given  the 
marginal  manner  in  which  the  MOGS  algorithm  performed,  given  all  of 
the  information  it  needed,  it  Is  doubtful  it  could  effectively  perform 
with  even  less  information. 

Beach  environments  of  the  type  being  studied  in  this  investiga¬ 
tion  (bare  soil,  no  vegetation)  are  typically  no  wider  than  one  to  two 
hundred  meters,  given  the  pixel  size  of  current  satellite  MSS  systems 
(approximately  80  meters),  it  is  unlikely  that  more  than  two  or  three 
pixel  wide  strip  of  information  could  be  obtained.  This  would  be  of 
very  little  use  in  beach  trafficabil ity  modeling.  The  improved  reso¬ 
lution  of  the  Landsat-D  and  SPOT  systems  offer  the  potential  to  better 
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Table  7.  Spectral  Parameters 


M0GS  Algorithm 

Landsat  1-3 

SPOT 

Landsat-D 

0.43-0.47 

0.47-0.49 

0.45-0.52 

0.49-0.51 

0.51-0.53 

0.56-0.59 

0.59-0.63 

0.5-0. 6 

.50-. 59 

0.52-0.60 

0.63-0.67 
0.70-0.75  • 

0.6-0. 7 

.61-. 69 

0.63-0.69 

0.75-0.80 

0. 7-0.8 

0.76-0.90 

0.80-0.90 

0.90-1 .00 
1.10-1.20 
1.20-1.35 
1.50-1.85 

0.8-1 .10 

.79-. 90 

1 .55-1 .75 

2.10-2.15 

2.08-2.35 

Spatial 

Resol ution 

ERIM  MSS8 
Scanner 

Landsat  1-3 

SPOT 

Landsat-D 

3  m  x  3  m 

57  m  x  79  m 

20  m 

30  m  x  30 

*Coastal  Zone  Color  Scanner 


CZCS* 

0.42-0.44 

0.50-0.5 

0.53-0.5 

0.65-0.6 

0.73-0.7 

1  .148-1  .152 


CZCS* 

825  m  x  825  m 
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discriminate  beach  environments,  while  the  coarse  resolution  of  the 
CZCS  system  offers  virtually  no  potential. 


Also  to  be  considered  in  using  satellite  systems  is  the  timeli¬ 
ness  of  the  data.  Currently,  data  from  a  specific  area  is  only 
available  once  every  several  days.  The  data  then  has  to  be  corrected 
for  atmospheric  effects.  In  summary,  it  is  felt  that  there  is  very 
little  potential  for  using  satellite  collected  data  in  beach  environ¬ 
ment  modeling  given  todays  operational  MSS  satellites. 
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APPLICATION  OF  A  MOGS  TYPE  ALGORITHM 
TO  NON-SAND  SOILS 

This  section  of  the  report  discusses  modeling  of  non-sand  soils. 
Specifically,  the  question  addressed  is  whether  a  MOGS  type  algorithm 
could  be  applied  to  mud  flats  where  the  water  to  particle  volume  ratio 
is  very  high. 

There  are  two  important  assumptions  required  by  the  reflectance 
model  of  moist  beach  sands  and  natural  water  bodies.  First,  one 
assumes  that  the  air-water  or  air-sand  boundary  is  essentially  planar. 
In  the  case  of  beach  sands,  this  planar  surface  is  obtained  naturally 
by  wind  on  dry  sand  or  by  rain  and  wave  action  on  wet  sand.  Gravity 
maintains  a  planar  surface  on  water  except  at  high  sea  states  when  the 
model  reflectance  predictions  will  begin  to  deviate  from  field 
measurements.  On  the  other  hand,  soils  containing  silt  and  clay  can 
sustain  a  non-planar  macrostructure  of  the  wet  or  dry  surface.  This 
macrostructure  will  tend  to  alter  the  directional  properties  of 
reflectance  by  significant  amounts.  Although  spectral  variations  will 
be  the  same,  the  magnitude  of  reflectance  in  a  given  direction  will 
change. 

Second,  soils  containing  silt  or  clay  tend  to  compact  with 
various  kinds  of  microstructures  depending  upon  soil  history.  A  basic 
model  assumption  concerning  the  scattering  of  radiation  by 
particulates  is  that  the  particulates  are  randomly  distributed 
throughout  each  stratum.  Optical  contact  between  particulates  is 
assumed  to  have  negligible  effect  in  the  scattering  process.  In  the 
case  of  sand,  the  microstructure  is  assumed  to  be  a  random  stacking 
with  the  area  of  points  of  contact  between  grains  negligibly  small 
compared  to  particulate  surface  area.  In  this  regard,  all  sands  have 
a  similar  microstructure  and  the  points  of  contact  do  not  influence 
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the  scattering  significantly.  In  saturated  sands,  the  scattering  is 
fundamentally  the  same  as  that  of  a  very  turbid  natural  water  body,  so 
turbid  with  sand  that  the  particulates  rest  upon  each  other  and  occupy 
a  major  fraction  of  the  volume  of  the  water  body.  As  water  is  drained 
to  field  capacity,  air  is  drawn  into  the  void  spaces  but  a  continuous 
film  of  water  surrounds  each  sand  grain.  The  model  treats  the  air  in 
the  voids  as  merely  a  new  kind  of  scattering  particulate  in  the  turbid 
water  with  a  total  volume  equal  to  the  total  volume  of  air  in  voids 
and  a  concentration  equal  to  the  number  of  voids  per  unit  sand  volume. 
The  number  and  size  of  voids  will  be  directly  related  to  the  number 
and  size  of  sand  grains  which  make  these  voids  by  random  packing. 

In  some  soils  containing  clay,  the  void  spaces  may  be  made  as  a 
result  of  preferential  platelet  to  platelet  stacking  induced  by  the 
surface  tension  of  the  water  adhering  to  the  particulates.  Therefore, 
the  size  of  voids  will  be  microstructure  dependent.  In  addition,  the 
areas  of  optical  contact  between  particulates  may  be  a  large  fraction 
of  the  area  of  a  particulate  and  thereby  favor  optical  transmission 
from  particulate  to  particulate  through  very  thin  water  interfaces 
instead  of  scattering  through  water  and  air  voids  to  arrive  at  the 
next  particulate  as  in  sands.  The  net  effect  will  be  a  scattering 
process  similar  to  that  for  a  much  larger  particulate  than  the 
Individual  clay  particles.  The  size  of  this  composite  particulate 
will  be  microstructure  dependent. 

The  proper  modeling  of  planar  surfaces  of  soils  would  require 
that  this  microstructure  be  taken  into  account.  Although  the  spectral 
features  would  be  predicted  correctly  by  the  present  sand  model,  the 
quantitative  accuracy  as  to  proportions  of  water,  air,  and  mineral 
would  not  likely  be  accurate  without  modification. 

The  modeling  of  non-planar  soil  surfaces  requires  the 
introduction  of  a  new  driving  parameter  which  describes,  in  some 
suitable  statistical  manner,  the  variations  of  the  surface  from  a 
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plane.  In  order  for  a  model  which  incorporates  surface  irregularity 
to  be  useful,  it  would  be  necessary  to  show  that  the  natural  soil 
surfaces  of  interest  become  irregular  in  a  relatively  few  character¬ 
istic  ways.  Otherwise,  there  would  be  an  unmanagably  large  number  of 

possibilities  for  directional  reflectance  variations  depending  upon 
chance. 


31 


APPLICATIONS  DIVISION 
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SUMMARY  AND  CONCLUSIONS 

The  four  year  activity  summarized  in  this  report  represents  a 
significant  contribution  to  the  understanding  of  the  spectral 
reflectance  response  of  beach  sands  in  the  .4  to  2.5  pm  region  of  the 
electromagnetic  spectrum.  The  spectral  shape  of  five  mineralogical 
types  were  studied  in  detail  and  a  radiative  theory  to  adequately 
describe  the  spectral  shapes  of  the  beach  sands  was  developed  through 
the  generation  of  the  "Suits  AQUASAND"  radiative  transfer  model.  This 
model,  an  adaptation  of  a  vegetation  reflectance  model,  calculates 
sand  reflectance  based  on  mineralogical  inputs  and  is  currently  being 
used  by  University  of  Michigan  faculty  and  scientists  from  the  United 
States  Geological  Survey. 


The  final  year  and  one-half  of  this  research  program  was 
concerned  with  development  of  a  methodology  whereby  an  algorithm  (the 
MOGS  algorithm)  was  produced  to  predict  mineralogy,  moisture,  and 
grain  size  of  unknown  beach  sands.  The  developed  MOGS  algorithm 
appeared  to  work  quite  well  in  prdicting  the  general  mineralogy  type 
and  moisture  of  sands  similar  to  the  ones  used  to  develop  the  moael , 
however,  difficulties  arose  when  trying  to  predict  grain  size  and 
classifying  sands  different  than  the  ones  used  to  produce  the  model. 
It  is  felt  that  the  shortcomings  of  the  MOGS  algorithm  could  be 
overcome  if  a  larger  number  of  beach  sand  samples  were  used.  The 
basic  methodology  used  to  develop  the  model  is  believed  to  be  sound 
and  can  be  used  in  other  remote  sensing  applications. 

One  application  for  the  MOGS  algorithm  was  to  determine  whether 
or  not  remote  sensing  data  can  be  used  in  beach  trafflcabil ity  models. 
Prior  to  the  development  of  the  MOGS  algorithm,  it  was  believed  that 
grain  size  and  soil  moisture  were  the  two  most  important  variables  in 
beach  trafficabil ity  models.  Recent  studies  (Thomson,  et  al . ,  1980) 
have  indicated  that  this  is  not  entirely  true. 


33 


I 


APPLICATIONS  DIVISION 


This  later  study  determined  that  of  all  the  factors  needed  to 
predict  a  soil  penetration  resistance  (cone  index)  for  a  given 
situation,  soil  moisture  and  jn,  situ  void  ratio  were  the  most 
important.  Grain  size  was  relatively  unimportant  in  the  sense  that 
only  gross  gradation  of  grain  sizes  are  needed  as  inputs  into  the 
beach  trafficabil ity  model.  It  was  also  reported  that  the  soil 
moisture  required  by  the  beach  trafficabil  ity  model  was  not  a  surface 
moisture,  as  determined  by  the  MOGS  algorithm,  but  a  soil  moisture 
profile  over  the  first  15  centimeters  of  the  sand.  The  MOGS  algorithm 
produces  a  soil  moisture  for  the  first  centimeter  of  sand  at  the 
surface. 

The  nain  conclusion  with  respect  to  beacn  trafficabil  ity  is  that 
in  orde"  to  use  present  day  beach  trafficabil ity  models,  additional 
information  is  needed:  a  third  dimension  (depth)  of  soil  moisture 
data  and  an  estimate  of  J_n  situ  void  ratio.  Up  to  the  present, 
neither  of  these  parameters  have  been  studied  in  detail  to  ascertain 
as  to  whether  or  not  they  can  be  remotely  sensed. 

Remote  sensing's  utility  in  providing  trafficabil  ity  inputs  is  in 
providing  synoptic  and  temporal  coverage  of  a  given  geographic 
location.  A  hybrid  approach  of  ground  truth  test  plots  coupled  with 
remote  sensing  data  to  determine  areal  extent  of  a  given  set  of 
environmental  parameters  is  perhaps  the  most  plausible  approach  in 
solving  the  trafficabil ity  input  problem. 

To  summarize  the  overall  four  year  research  activity,  it  can  be 
stated  that  remote  sensing  algorithms  using  passive  .4  to  2.5  um  data 
have  been  developed  to  predict  moisture  and  grain  size  for  a  limited 
class  of  beach  materials.  These  algorithms  can  provide  useful  infor¬ 
mation  for  a  variety  of  coastal  geographic  studies,  but  inherent 
limitations  of  the  technique  preclude  its  utilization  as  a  primary 
data  source  for  beach  trafficabil ity  predictions. 
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APPENDIX 


This  appendix  contains  three  publications  which  resulted  from  work 
done  for  this  contract.  These  papers  are: 

1.  "AQUASAND:  A  Beach  Reflectance  Model  and  Validation  Tests," 
published  in  the  Proceedings  of  the  Fifth  Canadian  Symposium  on 
Remote  Sensing  of  Environment,  Victoria,  B.C.,  1978. 

2.  "The  Use  of  Remote  Sensing  in  the  Determination  of  Beach  Sand 
Parameters,"  published  in  the  Proceedings  of  the  Thirteenth 
International  Symposium  on  Remote  Sensing  of  Environment, 

Ann  Arbor,  Michigan,  1979. 

3.  "Determination  of  Beach  Sand  Parameters  Using  Remotely  Sensed 
Aircraft  Reflectance  Data,"  accepted  for  publication  in  Remote 
Sensing  of  Environment,  1 981 . 
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AQUASAND:  A  BEACH  REFLECTANCE  MODEL 
AND  VALIDATION  TESTS* 

R.A.  Shuchoan ,  G.H.  Suits,  C.F.  Davis 
Environmental  Research  Institute  of  Michigan 
Ann  Arbor,  Michigan  4S107 


ABSTRACT 

A  new  sand  reflectance  model  called 
AQUASAND  has  recently  been  developed.  This 
model,  a  modification  of  the  Suits  radiative 
transfer  vegetation  canopy  directional  reflec¬ 
tance  modal,  accounts  for  the  reflectance  and 
transmittance  of  sand  parciculates  in  the  .35 
to  2.5  um  spectral  range.  The  AQUASAND  model 
also  accounts  for  Che  influence  of  soil  mois¬ 
ture  wichin  the  sand.  The  model  will  ulti¬ 
mately  determine  the  practical  limits  of  re¬ 
mote  sensing  algorithms  for  determining  physi¬ 
cal  and  chemical  properties  of  beaches. 

The  AQUASAND  model  uses,  as  inputs,  the 
coefficients  of  absorption  and  scattering  and 
the  forward  scattering  fraction  for  each  min¬ 
eral  comprising  the  u»ach  sand  (i.e.,  quarts, 
feldspar,  kaolinice,  etc.),  Che  average  number 
of  grains  (particles)  par  given  volume  from 
which  average  cross  sections  of  each  mineral 
type  can  be  computed,  void  space,  and  the 
moisture  depth  profile  of  the  beach  to  calcu¬ 
late  che  reflectance  of  the  beach.  ERIM, 
using  its  Cary  14  spectral  rnfleccometer,  mea¬ 
sured  the  hemispherical  transmittance  and  re¬ 
flectance  of  sixteen  coomon  beach  forming  min¬ 
erals  at  a  number  of  prescribed  thicknesses  to 
obtain  che  needed  coefficients  and  forward 
scattering  fraction. 

Additionally,  ERIM  used  che  same  instru¬ 
ment  to  measure  actual  sand  samples  collected 
from  fifty  diverse  beaches  found  on  United 
States  coastlines.  These  reflectance  spectra 
were  used  to  evaluate  che  AQUASAND  modal.  The 
results  are  very  encouraging,  showing  good 
agreement  between  che  actual  beach  spectra  and 
the  modal  results  when  Che  individual  mine-  d 
components  comprising  Che  beach  are  correctly 
identified  and  inputed  into  the  model. 

The  validated  AQUASAND  modal  (based  on  ten 
of  che  fifty  actual  beach  samples)  is  current¬ 
ly  being  used  to  generate  algorithms  that  pre¬ 
dict  grain  siaa,  moisture  contents,  and  min¬ 
eral  composition  using  remotely  sensed  reflec¬ 
tance  information  in  the  .35  to  2.5  ^m 
spectral  range. 


•This  work  la  supported  by  che  United  States 
Office  of  Naval  Research  (ONR),  Contract  No. 
N0014-74-0273.  Dr.  Hans  Dolezalek  Is  che 
technical  monitor. 


RESUME 

Un  nouveau  modile  du  facteur  de  reflexion  des 
plages,  dinommi  AQUASAND,  a  ice  mis  au  point 
rAcemmenc.  Ce  modile  qul  esc  une  variance  du 
modile  SUITS  du  facteur  de  reflexion  dlrecclf 
du  tranaferc  de  rayonnement  du  couverc 
vigitatlf,  prend  en  compce  les  facteurs  de 
reflexion  ec  de  transmlsalon  des  parclcules 
de  sable  dans  l'incervalle  spectral  de  0.35  & 
2.5  um.  Ce  modile  prend  egalement  en  consider¬ 
ation  l’eff-t  de  l'humiditi  dans  le  sable.  Le 
role  de  ce  modele  esc  de  determiner  en 
definitive  les  limltes  pratiques  des  algorith- 
mes  de  ceiedicection  servant  a  etablir  les 
proprietia  physiques  ec  chiolques  des  plages. 

Les  donnies  d'entrAe  introduces  dans  AOUASAND 
pour  le  ralcul  du  facteur  de  reflexion  des 
plages  sont:  les  coefficients  d'absorptlon  ec 
de  diffusion,  ainsi  que  la  fraction  de  pro- 
diffusion  de  chacun  des  minerals  entranc  dans 
la  composition  du  sable  de  plage  (quartz, 
feldspath,  kaolinice,  ecc.);  le  nombre  moyen 
de  particules  par  unite  de  volume,  donnAe  que 
l'on  utilise  pour  calculer  la  section  de  chaque 
type  de  mineral;  l'espace  des  vides;  ec  la 
repartition  en  profondeur  de  l'humldlte. 

ERIM,  au  moyen  de  son  rAf lectometre  spectral 
Cary  14,  a  meaure  le  facteur  de  transmission 
hemisphirlqua  ec  le  facteur  de  reflexion  de 
seize  des  princlpaux  minerals  entranc  dans  la 
composition  des  plages,  &  des  epaisseurs 
prescrltes  pour  determiner  les  coefficients 
indispensmbles  et  la  fraction  de  pro-diffusion. 

En  outre,  ERIM  s'est  servl  du  meme  instrument 
pour  etudler  dea  ichanclllons  de  sable  prAlevis 
dans  cinquante  plages  differences  des  llctoraux 
amSricains.  Lee  spectres  du  facteur  de 
rCflextlon  obtenua  one  Ate  employes  pour 
apprAcler  le  modile  AQUASAND.  Les  resulcats 
sont  tris  encourageancs ;  11  y  a  un  accord  cris 
bon  encre  les  spectres  da  plage  reels  ec  les 
resulcats  slmulAs  quand  les  minerals 
individuals  conadcuanc  la  plage  sont  blen 
identifies  et  Introduits  dans  le  modile. 

La  modile  AQUASAND  valldi  (reposant  sur  10  des 
50  echantlllons  de  plage  resile)  esc  utilise  v 
reguliirement  pour  Alaborer  des  algorlchmes 
qul  prAdisent  la  granulomAtrle,  les  censure 
en  humlditi  ec  la  composition  en  minerals,  en 
se  servant  du  facteur  ds  reflexion  teiidAcecte 
dans  l'lntervalle  du  spectra  comprls  encre  0.35 
et  2.5  um. 

1.0  INTRODUCTION 

This  paper  discusses  a  new  radiative 
transfer  model  called  AQUASAND.  This  model 
predicts  the  reflectance  of  beech  sands  in  the 
.35-2.5  ftu  wavelength  range  of  the  electromag¬ 
netic  spectrum.  This  newly  developed  sand 
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reflectance  modal  will  ulcisutely  determine 
che  practical  limits  of  remote  sensing  algo¬ 
rithms  for  dataraining  physical  and  chaaical 
properties  of  beaches. 

The  use  of  .35  to  2.5  um  radiation  to 
determine  physical  and  chaaical  properties  of 
beaches  is  deeaed  feasible  froa  che  proven 
capability  of  mapping  silicate  minerals  using 
rescscrahlen  techniques  in  che  3  to  14  ua 
spectral  range,  a  procedure  developed  by 
Vincient  and  Thomson  (1972)  and  earlier  by 
Vincienc  and  Hunt  (1958).  The  .35  to  2.5  ua 
region  of  Che  speccrua  is  of  interest  for 
beach  paraaMCer  senaing  since  water  is 
relatively  transparent  in  this  region  (unlike 
che  8  to  14  ua  region)  and  beaches  typically 
have  high  moisture  concents. 

2.0  MODEL IHC 

The  capability  of  determining  important 
physical  properties  of  beach  sanda  by  reacts 
senaing  techniques  depends  upon  che  interac¬ 
tion  of  radiation  with  the  conatituent  materi¬ 
als  of  th*  beach.  Such  interaction  ia  cer¬ 
tainly  coaplex  but,  nevertheless,  must  follow 
the  laws  of  nature.  The  purpose  of  making  a 
mathematical  reflectance  model  of  Chia  coaplex 
interaction  is  to  achieve  insight  into  Che 
relationship  between  Che  remotely  received 
signals  and  che  physical  properties  of  the 
beach  sand  chac  are  of  interest. 

A  mathematical  model  of  a  physical  phenom¬ 
enon  is  che  result  of  incorporating  che  mathe¬ 
matical  expressions  of  che  laws  of  nature  aa 
they  apply  to  a  coaplex  circuautance  so  that 
che  conclusions  drawn  froa  che  assembled  ex¬ 
pression  correspond  to  the  physical  results 
of  an  experiment  under  similar  circumstances. 
Such  a  modal  requires  experimental  validation 
using  a  few  but  diverse  circumstances  in  order 
to  prove  chac  Che  essence  of  the  phenomenon  is 
contained  in  che  logic  of  che  modal.  Once 
validation  has  been  achieved,  che  aodel  may  be 
used  as  a  readily  available  and  inexpensive 
substitute  for  experiments  under  all  circum¬ 
stances  within  Che  scope  of  Che  circumstances 
of  che  validation  experiments.  In  addition, 
che  logical  structure  of  che  model  provides 
che  insight  into  che  significant  interaction 
processes  so  chac  more  general  conclusions  aay 
be  drawn. 

2.1  Reflectance  Model  Concepts 

The  most  elementary  aodel  of  sand  reflec¬ 
tance  is  the  simple  plana  mixtures  aodel.  The 
aodel  employs  the  assumptions  that  all  sand 
particles  are  opaque  and  are  randomly  mixed. 
The  surface  of  the  sand  layer  exposes  sand 
particles  in  proportion  to  the  product  of 
their  mean  cross  sectional  area  and  the  con¬ 
centration  of  the  particles  in  the  sand  mix¬ 
ture.  Thus,  the  reflectance  spectrum  of  the 


mixture  is  predicted  to  be  the  area  weighted 
average  of  the  reflectance  spectra  of  che  var¬ 
ious  minerals  exposed  at  the  sand  surface. 
Multiple  scattering  between  particles  is 
assumed  to  be  negligible  and  one  surface  par¬ 
ticle  ia  asausMd  not  to  obscure  froa  view  an 
adjacent  surface  particle. 

This  eleaentary  aodel  fails  to  achieve 
good  accuracy  because  the  transmittance  of 
particles  in  a  finely  divided  state  may  be 
quite  large  and  multiple  scattering  of  radia¬ 
tion  between  particles  should  be  significant. 
In  addition,  the  packing  of  grains  in  layers 
produces  partial  exposure  at  che  surface  so 
chat  line  of  sight  co  some  particles  will  de¬ 
pend  upon  Che  direction  of  view.  A  more  com¬ 
plex  model  is  required  to  account  for  these 
effects. 

The  more  complex  model  which  is  used  in 
this  work  employs  the  identical  concepts  that 
are  employed  by  the  directional  reflectance 
model  for  vegetative  canopies  (Suits,  1972). 

It  may  be  hard  co  visualise  off-hand  that  the 
interaction  of  radiation  with  a  vegetative 
canopy  is  homomorphic  with  Che  interaction  of 
radiation  with  sands  because  vegetation  and 
sand  hardly  appear  the  same  to  the  eye. 
Nevertheless,  Che  essence  of  the  reflection, 
transmission,  and  multiple  scattering  phenom¬ 
ena  ia  the  same  concept  for  the  optical  param¬ 
eters  of  the  components  that  are  involved  aa 
long  aa  the  wavelength  is  much  smaller  than 
che  particles. 

The  first  assumption  of  this  model  is  chac 
che  scattering  components  are  distributed  son 
or  less  uniformly  in  horizontal  layers  where 
Che  mixture  of  component  types  may  be  differ¬ 
ent  in  the  various  layers.  In  a  corn  field, 
for  example,  Che  tassels  slwaya  appear  at  Che 
cop,  green  healthy  mature  leavea  appear  in  a 
middle  layer  and  necrotic  leaves  appear  usual¬ 
ly  near  Che  soil.  In  sands,  the  action  of 
wind  and  waves  Is  likely  to  stracify  mineral 
mixtures  vertically  and  certainly  moisture 
content  varies  with  a  vertical  moisture  pro¬ 
file.  The  division  in  layers  is  done  in  order 
to  quantify  statistically  what  may  otherwise 
be  nearly  a  continuous  distribution. 

The  second  assumption  is  chat  che  radia¬ 
tion  field  may  be  divided  into  two  types  of 
radiant  flux  —  specular  and  diffuse.  The 
specular  flux  represents  the  radiation  arriv¬ 
ing  from  Che  source  with  rectilinear  propaga¬ 
tion  and  passes  through  che  holes,  cracks,  and 
voids  of  the  ensemble  of  randomly  packed  com¬ 
ponents  without  deviation.  The  diffuse  flux 
is  derived  from  Che  specular  flux  and  ia  that 
part  of  the  specular  flux  which  has  been  in¬ 
tercepted  by  a  scattering  component  ac  least 
once  and  ia  scattered  in  both  forward  and 
backward  directions.  In  a  vegetative  ces-py, 
specular  flux  frequently  reaches  the  soil 
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Level  end  appears  it  sun  flecks  on  Che  soil. 

In  e  send,  speculer  flux  diminishes  exponen- 
tislly  Co  negligible  proporcions  in  only  e  few 
millimeCers  depch.  Opcictlly,  Che  send  is 
infinicely  deep.  However,  Che  diffuse  flux  is 
derived  from  Che  speculer  flux  in  Che  idcnti- 
cel  meaner.  The  diffuse  flux  may  peneCrsCe 
much  deeper  inco  Che  send  chen  cen  speculer 
flux. 

The  chird  model  essuapcion  is  chsc  Che 
menner  of  scsccering  by  minerel  perciculeees 
csn  be  adequately  repretenced  by  replecing 
etch  minerel  perciculece  wich  a  sec  of  equiva- 
lenc  Lembercien  penels  which  heve  Che  seme 
speccrel  crentmiccence  end  refleccence  es  does 
che  componenc.  This  assumption  defines  a  sim¬ 
plified  scaCCering  phase  function  which  per¬ 
mits  one  Co  calculate  ensemble  reflectances  in 
closed  form.  The  form  of  che  scattering  phase 
function  becomes  significant  when  tingle  scat- 
cering  is  che  dominant  phenomenon.  Scattering 
by  widely  dispersed  aerosols  in  che  atmosphere 
(e.g.,  smoke  and  dust)  exhibits  detailed  phase 
function  effects.  However,  as  che  degree  of 
multiple  scattering  increases,  che  detailed 
feacurea  of  che  scattering  phase  function  are 
no  longer  significant.  In  both  sands  and 
vegetative  canopies,  multiple  scattering 
effects  dominate  because  of  the  high  density 
of  scattering  components. 

The  fourth  model  assumption  is  thac  the 
diffuse  flux  moves  generally  vertically  upward 
and  downward  wich  a  Lambertian  angular  distri¬ 
bution  as  a  first  approximation.  The  reflec¬ 
tance  of  Che  ensemble  is  calculated  using  che 
method  of  self-consistenc  field  wich  che  spec¬ 
ular  and  approximate  diffuse  flux  as  che  illu- 
minant  of  components.  The  ensemble  reflec¬ 
tance  is  not  necessarily  Lambertian  but  is 
only  approximately  so.  Both  vegetation  and 
sand  meet  this  approximate  criterion. 

Because  of  che  homomorphic  relationship  of 
radiative  interactions  with  sand  and  vegeta¬ 
tive  canopies,  che  reflectance  modal  previous¬ 
ly  developed  for  vegetation  is  applied  to  sand 
wich  che  appropriate  component  properties  for 
sand  minerals  substituted  for  vegetative 
canopies. 

2.2  Optical  Properties  of  Sand  Components 

The  cross  section  view  of  a  hypothetical 
sand  layer  containing  two  kinds  of  minerals  is 
shown  in  Figure  1.  The  irregular  shapes  of 
che  sand  grains  result  in  some  store  or  lass 
random,  loose  packing  with  many  voids.  These 
grains  are  to  be  replaced  by  a  number  of 
equivalent  Lambertian  panels  which  will  inter¬ 
cept  approximately  the  same  amount  of  radiant 
flux  as  do  che  actual  grains.  The  spectral 
reflectance  and  transmittance  of  Che  panels 
are  to  be  the  seme  as  Che  spectral  reflectance 


and  Cransmittance  of  Che  grains.  Since  Chese 
speccral  properties  may  change  with  grain  size 
or  state  of  division,  some  means  of 
calculating  Che  appropriate  speccral  proper- 
cies  is  required. 

An  auxiliary  refleccance  and  cransmictance 
model  for  mineral  Chin  sections  was  developed 
in  order  co  relate  the  inherent  spectral  prop¬ 
erties  which  are  characteristic  of  a  mineral 
co  che  properties  of  chat  mineral  in  any  state 
of  division.  Three  characteristic  bulk  prop¬ 
erties  were  caken  to  be  sufficient  for  this 
purpose  —  Che  speccral  absorpeion  coeffi¬ 
cient,  a,  Che  forward  scattering  fraction,  FS, 
and  che  scaCCering  coefficient,  s. 

Radiation  which  penetrates  a  mineral  may 
be  absorbed  and  converted  into  heat  energy, 
depending  upon  che  chemical  composition  of  Che 
mineral.  A  mineral  which  is  inCemally  homo¬ 
geneous  wichout  inclusions  and  cracks  will 
transmit  radiation  passing  chrough  it  in 
accordance  wich  the  relation, 

ECx)  -  E0  e*«,  (1) 

where  E(x)  is  Che  irradiance  on  a  plane  aC 
depch  x  in  the  mineral, 

E0  is  che  irradiance  on  a  plane 
inside  the  first  surface, 

a  is  Che  spectral  absorption 
coefficient. 

The  spectral  absorption  coefficient  will 
be  a  function  of  Che  wavelength  of  che  pene¬ 
trating  radiation  and  will  depend  upon  che 
chemical  composition  of  che  mineral.  The 
speccral  absorpeion  coefficient  is  largely 
responsible  for  the  speccral  variations  in 
mineral  refleccance  and  cransmittance. 

The  formation  of  minerals  is  a  complex 
natural  process  so  thac  minerals  may  not  be 
optically  homogeneous.  Foreign  materials  are 
often  formed  in  che  interior  of  che  mineral. 
Such  features  as  fractures,  gas  bubbles,  small 
crystals  of  associated  minerals,  and  grain 
boundaries  of  anisotropic  crystals  create  in- 
homogenaicies  within  mineral  bodies.  These 
inhomogeneities  reflect  or  scatter  and  deviate 
penetrating  radiation  from  rectilinear  propa¬ 
gation.  A  collimated  beam  of  radiacion  which 
propagates  rectilinear ly  chrough  che  body  of  a 
mineral  will  be  diminished  due  co  such  scaC¬ 
Cering  by  che  relation, 

E(x)  -  E0  a-**,  (2) 

where  E(x)  is  che  irradiance  of  rectilinear 
flux  on  a  plane  at  depch  x, 

E0  is  the  irradiance  of  rectilinear 
flux  at  the  first  surface, 

s  is  che  scattering  coefficient. 
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The  radiaeion  which  is  scattered  will  tend 
to  propagate  either  deeper  into  the  mineral  — 
forward  scatter  —  and  contribute  to  transmit¬ 
tance  or  reverse  direction  and  propagate  back 
out  of  the  mineral  —  backacatter  —  and 
contribute  to  reflectance  of  the  mineral.  The 
fraction  of  scattered  radiation  which  contin¬ 
ues  deeper  into  the  mineral  is  the  forward 
scattering  fraction(  FS.  The  fraction  of 
scattered  radiation  which  reveraaa  ia  then 
(1-FS).  The  scattering  coefficient  and  the 
forward  scattering  fraction  can  ba  spectrally 
dependent  but  should  not  be  the  primary  deter¬ 
minant  of  the  spectral  quality  of  a  mineral. 

These  three  optical  properties  of  a  min¬ 
eral,  a,  s,  and  FS,  are  assumed  to  be  inde¬ 
pendent  of  che  thickness  of  the  mineral.  That 
is,  the  inhomogeneous  structure  and  chemical 
composition  of  a  mineral  thin  section  is 
assumed  to  be  evenly  distributed  so  Chat  a,  s, 
and  FS  are  inherent  properties  characteristic 
of  the  kind  of  mineral  and  not  che  size  of  Che 
mineral  sample. 

The  auxiliary  reflectance  and  transmit¬ 
tance  model  for  mineral  chin  sections  makes 
use  of  chese  three  properties  and  the  index  of 
refraction  of  che  mineral  to  yield  che  thin 
section  transmittance  and  reflectance  for  any 
mineral  thickness.  The  value  of  a,  s,  and  FS 
for  each  wavelength  must  be  determined  experi¬ 
mentally.  However,  these  properties  cannot  be 
determined  by  direct  experiment.  Instead,  a 
spectral  reflectometer  was  used  to  measure  the 
transmittance  and  reflectance  of  mineral  chin 
sections  having  various  thicknesses.  Since  a, 
s,  and  FS  are  presumably  thickness  invariant, 
che  value  of  a,  s,  and  FS  may  be  determined  by 
finding  che  value  of  a,  s,  and  FS  which,  when 
used  in  calculating  che  cranamictances  and 
reflectances  of  thin  sections  having  these 
various  thicknesses,  yield  matching  values  for 
che  reflectances  and  transmittances  found  ex¬ 
perimentally.  A  computer  iteration  technique 
was  used  for  this  purpose. 

The  spectral  values  of  a,  s,  and  FS  for 
some  of  the  consnon  sand  minerals  have  been 
tabulated.  For  certain  opaque  mnerals,  such 
aa  hematite  and  limonite,  Che  valuea  of  a,  s, 
and  FS  could  not  be  determined.  The  transmit¬ 
tance  and  reflectance  of  these  minerals  are 
independent  of  grain  thickness  for  all  thick¬ 
nesses  that  are  likely  to  be  found  in  sanda. 
The  reflectance  and  transmittance  of  pure, 
unstained  clear  quartz  ia  due  almost  entirely 
to  surface  effects  which  are  also  independent 
of  grain  thickness. 

The  optical  propertiea  of  sand  components 
are  determined  and  are  introduced  into  che 
sand  reflectance  model.  The  spectral  trans¬ 
mittance  and  reflectance  of  the  equivalent 
Lambertian  panels  for  each  mineral  are  deter¬ 


mined  using  the  auxiliary  reflectance  and 
transmittance  model  for  mineral  thin  sections 
where  the  thickness  of  thin  section  is  Che 
mean  grain  thickness  for  each  mineral.  The 
mean  cross  section  of  grains  of  each  mineral 
type  is  multiplied  by  the  corresponding  number 
of  such  grains  per  unit  volume  for  a  given 
beach  sample  and  represents  the  scattering 
effect  of  the  equivalent  Lambertian  panels. 

2.3  Model  Description 

Figure  2  indicates  the  flow  of  information 
beginning  with  the  raw  experimental  CARY 
transmittance  (r)  and  the  reflectance  (*>)  data 
and  ending  with  the  predicted  reflectance 
spectrum  for  sand.  As  indicated  previously, 
the  needed  inputs  into  the  AQUASAND  model  are 
transmittance,  r,  and  reflectance,  p,  values 
for  the  individual  minerals  in  a  particular 
sand  configuration.  These  r  and  p  values  must 
be  calculated  for  each  particle  size.  Exp-'  - 
mentally,  r  and  p  spectra  were  determined 
each  miheral  at  given  thicknesses.  Some  min¬ 
erals  yhich  were  opaque  exhibited  a  negligible 
asiount  of  transmittance.  Iron  stains  were 
treatea  as  separate  minerals  and  experimen¬ 
tal  r  and  p  spectra  were  determined  for  them. 
These  inputs  are  shown  on  Figure  2. 

In  order  to  predict  r'sandp's  for  any  par¬ 
ticle  size,  properties  independent  of  size  can 
determined.  These  properties,  previously  dis¬ 
cussed,  are  the  forward  scattering  (FS),  scat¬ 
tering  coefficient  (a)  and  the  spectral  ab¬ 
sorption  coefficient  (a).  The  computer  pro¬ 
gram  calculates  these  coefficients  for  a  given 
mineral  using  the  experimental  r's  ar.j  p's  and 
the  experimental  thicknesses  and  the  param¬ 
eters  of  refraction  (particle  to  air  (RPA)  and 
air  to  particle  (RAF)). 

These  determinations  of  the  fundamental  a,  ' 
FS,  and  s  properties  of  sixteen  coaon  beach 
forming  minerals  in  the  .35-2.3  urn  range  are 
the  first  such  measurements  made  to  date. 

Figure  3  is  a  plot  of  a  (a)  and  s  for  a  car¬ 
bonate  beach. 

Once  the  absorption  and  scattering  coeffi¬ 
cients  are  determined,  another  computer  pro¬ 
gram  uses  these  coefficients  coupled  with  the 
parameters  RFA  and  RAF  to  predict  r's  and  p's 
for  a  given  thickness.  These  predicted  values 
are  then  used  to  model  the  sand. 

The  final  program,  AQUASAND,  computes  the 
reflectance  of  mixtures  of  minerels.  The  in¬ 
puts  to  this  program  consist  of  the  predicted 
r’sandp's  for  the  minerals  constituting  the 
sand,  the  amount  of  iron  stain,  the  depth  of 
the  sand,  the  reflectance  of  an  infinite  depth 
beckground  underneath  the  sand,  the  particle 
cross  sections  and  packing  factors  for  che 
various  minerals,  the  angle  of  illumination 
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from  a  point  source  illuminator  (or  optional 
Lambertian  iourca)  and  tha  angular  poaition  of 
tha  aanaor  and  tha  percent  noiatura  of  tha 
tand  on  a  millimeter  acala. 

Tha  output  of  AQUASAND  ia  a  pradietad  ra- 
flactanca  apactra  of  tha  aand  for  varioua  con- 
figurationa  of  illuminator  poaition  (0)  and 
aanaor  poaition  (0,  f ) .  Tha  AQUA SAND  program 
wee  darivad  from  a  program  which  pradietad 
raflactancaa  for  a  two  layered  vagatation  can¬ 
opy  with  a  background.  Thua  tha  aand  can  hava 
two  layara  with  diffarant  danaitiaa  and  grain 
aizaa  in  aach  layar.  A  background  macarial 
for  tha  aand  can  alao  ba  apacifiad.  Thaaa 
propartiaa  of  AQUASAND  allow  graat  flexibility 
for  modeling  varioua  configuration  auch  ar 
thin  aanda,  aand  after  a  atom  with  a  layar  of 
organic  utter  or  atonaa  on  top,  ate. 

3.0  RESULTS 

Sevan  beachaa  of  divaraa  mineralogy,  moia- 
ture  content,  and  grain  aiza  ware  aalactad  for 
evaluation  of  the  AQUASAND  nodal.  Tha  beachaa 
aalactad  for  model  verification  along  with 
their  mineralogical  and  phyaicnl  parametera 
uaad  aa  inputed  into  AQUA  SAND  are  givan  in 
Tab  la  1.  Figure  4  ahowa  the  aevan  AQUA  SAND 
generated  reflectance  curvea  (in  the  .4-2.3  urn 
region  of  the  apectrum)  for  tha  reapectiva 
beachaa.  Super inpoaad  on  aach  of  tha  aavan 
curvea  appeara  the  actual  apectrum  of  the 
beach  aa  aaaaurad  uaing  the  ERIM  Cary  14 
apactroraf lactometer. 

Thaaa  raaulta  are  vary  encouraging,  for 
although  the  abaoluta  valuea  of  tha  curvea 
differ,  tha  overall  ahape  and  abaorption  banda 
a  a  am  to  correapond  remarkably  wall. 

4.o  sintuurc 

The  AQUA  SAND  reflectance  modal  ha  a  bean 
aaciafactorily  evaluated.  Tha  andel  ia  cur¬ 
rently  being  uaad  to  generate  hypothetical 
aanda,  ao  that  grain  aiza  and  moiatura  predic¬ 
tion  algorithma  can  be  generated  uaing  remote 
aanaad  data. 

A  beach  reaenbling  tha  aineralogy  of  beach 
aanple  AD2  (aaa  Table  1,  and  Figure  4)  haa 
bean  parametrically  run  with  varioua  moiatura 
and  grain  aizaa  uaing  AQUASAND.  Tha  generated 
apactral  grapha  aimilar  to  thoaa  ah own  in 
Figure  4  ware  then  examined  to  determine  opti¬ 
ma  tingle  or  ratio  apactral  intervale  to  pre¬ 
dict  moiatura  and  grain  aizp. 

A  aumaary  of  that  analyeia  ia  ahowa  in 
Figure  3.  The  graph  on  the  left  of  the  figure 
ahowa  ratio  channela  aelected  while  the  grapha 
on  the  right  ahow  tingle  channela  aelected 
from  AQUASAND  parametrically  varied  raaulta. 


Currently  AQUASAND  ia  being  uaad  to  gener¬ 
ate  beach  aanda  where  only  one  phyaical  or 
chemical  property  ia  parametrically  varied. 
Ragreaaion  analyeia  ia  then  applied  to  the 
AQUASAND  raaulta  to  obtain  tha  finalized  moia- 
Cure  and  grain  aiza  prediction  algorithma. 
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Cross  Section  of  Send  Layers  EquiUrant  Lambertian  Psnsls 

Layer  1  is  mown  as  consisting  of  pains  of  only  one  kind  of  mineral. 

Layer  2  Is  shoem  to  consist  of  a  mixture  of  two  kinds  at  minerals. 

The  equivalent  Lambertian  scattering  panels  are  illustrated  on  the  right. 

FIGURE  1.  CROSS  SECTION  OF  SAND  LAYERS 


TABLE  1.  BEACH  SAND  CONSTITUENTS 


3«°d  Mean  Crain  Send 


l .  Carbonate 

99Z  Carbons ce 

11  Organics 

1  ™ 

2.  D-Dune  HI 

95Z  quarts 

31  Feldspar 

.19 

3.  MX2 

98Z  Quarts 

2Z  Carbonate 

.22 

4.  BA1 

98Z  Quarts 

2Z  Feldspar 

.40 

3.  EA2 

90Z  Quarts 

6Z  Feldspar 

2Z  Kao  Unite 

2Z  Opaques 

.28 

6.  AD2 

98Z  Quarts 

2Z  Feldspar 

.38 

7.  HA1 

35Z  Quarts 

28Z  Feldspar 

4Z  Kao  Unite 

13Z  Opaques 

.32 
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FIGURE  2.  DIAGRAM  OF  MODEL  FLOW 
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REFLECTANCE  ALGORITHM 


EVALUATION  OF  HICH  RESOLUTION  SIDE  LOOKING 
AIRBORNE  RADAR  ON  THE  UNIVERSITY  OF  CUELPH 
TEST  STRIP 

Srian  Brisco  and  Richard  ?rocz 
Department  or  Land  Resource  Science 
University  or  Guelph 
Guelph,  Ontario.  NIG  2W1 


ABSTRACT 

Imagerv  obtained  from  che  four  channel  ERIM 
radar  system  was  evaluated  for  soil  survey 
and  field  type  identification  purposes.  The 
area  imaged  was  snow  covered  to  a  maximum 
depth  of  approximately  1  meter.  An  icv  snow 
surface,  high  moisture  concent,  or  ice  layers 
wichln  che  snovpack  may  have  been  che  cauae  of 
attenuation  of  the  X-HH  signal,  creating  a 
lack  of  tonal  variations.  The  ocher  three 
channels  exhibited  tonality  patterns  which 
varied  between  and  within  channels  and 
agricultural  fields.  These  conalicy  patterns 
were  not  relaced  to  geologic  or  topographic 
Influences,  but  low  correlations  with  Che 
soils  map  of  the  area  were  observed. 

Vegetation  at  these  locations  appeared  co  be 
the  cause  of  these  correlations.  Grey  scale 
values  of  che  seven  dominant  field  types 
encountered  in  che  study  area  Indicated  chat 
discrimination  la  possible  using  tonal 
signatures  from  all  four  channels  of  this  SAR 
system.  The  cross-comparison  of  each 
channel's  response  value  may  separate  fields 
with  very  similar  averaged  reflectance  values. 
More  ground  truth  and  more  effective  use  of 
ancillary  data  are  needed  to  improve  che 
Interpretation  procedure.  A  more  quantitative 
method  of  che  tonal  measurement  is  also 
needed  and  combined  with  a  measure  of  texture 
would  greatly  enhance  che  discriminative 
ability  of  radar  systems.  The  combination  of 
these  improvements  would  make  possible  che 
digital  classification  of  SAR  data  by 
computer.  Thus  multi- temporal  and  multi- 
speccral  radar  imagery,  in  combination  with 
the  effective  use  of  ancillary  data,  have 
che  potential  co  discriminate  crops  and  soils 
for  agricultural  purposes. 


RESUME 

Des  Images  de  terrain  enneigi,  obcenuea  avec  un 
svstime  radar  de  quaere  bandes,  ont  ici 
ivaluies  en  vues  d'utillsaclon  pour  l'inventaire 
des  sols  ainsl  que  pour  la  reconnaissance  des 
caractirlstlques  des  surfaces  agricoles.  Une 
surface  glacie,  une  teneur  en  eau  ilev<e,  ou 
des  couches  de  glace  1  l'lneirleur  de  la 
neige  one  pu  causer  1' attinuacion  du  signal 


Presented  to  the  5th  Canadian  Symposium  on 
Remote  Sensing,  Victoria,  August  1978. 


X-HH  creant  ainsl  un  manque  de  variations  dans 
la  conalite.  Dans  les  crols  autres  canaux  les 
tonalieis  des  couleurs  varlaient  encre  les 
bandes  et  suivanc  le  type  de  champs  agricole. 

La  conaliti  des  couleurs  n'apparaic  pas 
ecre  fonction  de  facteurs  giolcglques  ou 
topographiques  ec  les  correlations  avec  le 
relev*  des  sols  pour  la  region  sonc  pauvres. 
Cecl  pourrait  ecre  du  aux  Influences  du  couvert 
vegetal.  Les  valeurs  obcenues  dans  l'echelle 
des  grls  pour  les  sept  types  de  champs 
indiquenc  qu'une  distinction  esc  possible 
suivanc  les  caracterlstiques  conalos  de  chacune 
des  quaeres  bandes.  Er  comparanc  les  valeurs 
obcenues  pour  chaque  bands,  des  champs  avec 
dec  reflectances  "moyennes"  tres  semblable 
peuvenc  etre  separfs  .  Ce  procedi  de  discrim¬ 
ination  pourrait  acre  amellore  par  une  mesure 
plus  quantitative  du  con  ec  une  meilleure 
reconnaissance  sur  le  terrain.  Ce  procede 
pourrait  amener  une  classification  automaclque 
par  ordlnaceur  des  donnees  obcenues  avec  un 
svseeme  SAR.  Une  telle  lmagerie  au  radar 
multi-cemporelle  ec  mulcl-speccrale,  utilisee 
effecclvement  avec  des  donnies  auxllialres 
peut  permettre  d' identifier  les  recoltes  et 
les  sols  agricoles. 

INTRODUCTION 

Tonal  variations  observed  on  radar  Images 
of  the  earth's  surface  depict  che  ability 
of  the  ground  targets  co  reflect  electro¬ 
magnetic  radiation  incident  on  che  Imaged 
targets.  The  proportion  of  incident  radi¬ 
ation  returned  co  the  radar  sensor,  i.e., 
backscactsr,  depends  on  dielectrical  and 
geometrical  properties  of  the  ground 
elements.  Various  research  efforts  in  the 
past  have  attributed  che  cone  differences 
observed  on  radar  imagery  co  characteristics 
of  snow  cover  (Waite  and  MacDonald,  1970; 
Vickers  and  Rosa,  1972),  surface  soil  condi¬ 
tions  (Morain  and  Campbell,  1974;  Cihlar  et 
aL,  1975),  and  crop  type  (Haralick  et  al. , 

1970;  Bush  and  Ulaby,  1975).  Radar  is  a 
sensor  of  great  potential  co  many  kinds  of 
studies  in  agronomy,  but  more  research  is 
needed  on  che  interactions  of  che  earth's 
surface  and  che  radar  waves  co  develop 
effective  Interpretation  procedures. 

With  che  advent  of  synthetic  aperture  radar 
systems  (SAR),  finer  resolution  Imagery  has 
provided  more  detailed  representations  of 
the  earth's  surface.  Dual  polarization, 
multi-frequency  radar  systems  are  che  micro- 
wave  equivalents  of  mulclspectral  scanners, 
which  have  proven  useful  for  many  remote 
sensing  applications  (Haralick  ec  al.,  1970). 
The  purpose  of  this  study  was  to  evaluate 
four  channel  SAR  Imagery,  of  snow  covered 
terrain,  for  information  pertaining  co  soil 
survey  and  cover  type  discrimination. 
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THE  USE  OF  REMOTE  SENSING  IN  THE  DETERMINATION 
OF  BEACH  SAND  PARAMETERS* 


C.F.  Davis,  R.A.  Shuchman,  and  G.H.  Suits 

Environmental  Research  Institute  of  Michigan 
Ann  Arbor,  Michigan 


ABSTRACT 

A  reflectance  model  (AQUASAND)  was  developed  to  gain  insight 
into  the  spectral  effects  of  changing  mineralogy,  moisture,  and 
grain  size  related  to  beach  sands.  The  model,  a  modification  of 
the  Suits  radlacive  transfer  vegetation  canopy  model,  uses  the 
transmittance  and  reflectance  of  the  component  minerals,  the  desired 
sand  moisture  content,  and  the  desired  sand  grain  size  to  produce  a 
reflectance  spectrum  in  the  .35  to  2.5  urn  range. 

Using  information  from  AQUASAND,  a  mineralogy,  moisture,  and 
grain  size  predictive  algorithm  (MOOS)  was  developed  based  on 
laboratory  sand  reflectance  spectra  Input  to  the  .algorithm  is  a 
set  of  17  spectral  reflectance  "bands"  simulating  an  airborne 
multispectral  scanner  (MSS)  configuration.  The  MOGS  algorithm  first 
determines  the  appropriate  mineralogical  class  for  the  input  sand. 
Within  a  given  class,  regression  equations  are  used  to  determine 
moisture  and  grain  size. 

The  predictive  results  of  the  MOGS  algorithm  are  very  encouraging. 
When  tested  on  70  of  the  sand  reflectance  spectra  from  which  it  was 
derived,  the  correlation  of  actual  to  predicted  moisture  was  967.. 

The  correlation  of  actual  to  predicted  grain  size  based  on  46  samples 
was  8u7..  Tests  on  independently  collected  sand  spectra  yielded  similar 
results.  The  algorithm  was  also  successfully  applied  to  actual  MSS 
data,  collected  over  the  Lake  Michigan  coastline,  to  generate  moisture 
distribution  and  grain  size  distribution  digital  images  of  the  beach 
region. 


1.  INTRODUCTION 

Historically  remote  sensing  has  been  proven  useful  in  the  delineation  of 
rock  types  most  often  with  the  applications  of  future  mineral  exploration  in 
mind.  In  this  investigation,  beach  sands  were  analyzed  with  the  Intention  of 
determining  not  only  mineralogy  but  also  moisture  and  grain  size.  These  three 
parameters  are  of  Interest  from  a  beach  traff icability  and  sediment  transport 
point  of  view.  Recognition  of  these  parameters  would  also  allow  the  identifica¬ 
tion  of  beach  mineral  deposits  based  on  grain  size  and  mineralogy.  Since 
beaches  are  formed  by  inttnse  erosion  of  the  parent  materials,  harder  minerals 
tend  to  be  preserved  and  concentrated  while  others  are  dissipated.  Depending 
on  the  origin  of  the  sand,  these  residual  minerals  may  be  economically  valuable. 

In  order  to  fully  understand  the  effects  of  changes  in  mineralogy,  moisture, 
and  grain  size  on  the  spectral  reflectance  of  beach  sand,  a  sand  reflectance 
model  was  used.  The  model,  an  adaptation  of  the  Suits  radiative  transfer 
vegetation  canopy  model,  is  known  as  AQUASAND  (Suits,  1972).  It  uses,  as  inputs, 

^This  work  is  supported  by  the  Office  of  Naval  Research  (ONR)  Contract  No. 

N0014-74-C-0273.  Mr.  Hans  Doiezalck  is  the  ONR  technical  monitor. 
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the  reflectance  and  cransmictance  for  each  mineral  comprising  the  beach  sand 
i.e.,  quartz,  feLdspar,  kaolinite,  etc.).  In  addition,  input  to  the  model 
• ncludes  the  sand  grain  size,  the  void  space  in  the  sand,  and  the  moisture  pro¬ 
file  as  a  function  of  depth.  3y  varying  these  inpuc  parameters  insight  was 
g  lined  into  the  effect  of  physical  changes  on  the  bulk  sand  reflectance. 

Using  the  information  obcainad  from  AQUASAND,  the  mineralogy,  moisture,  and 
grain  size  (MOCS)  algorithm  was  developed  using  reflectance  spectra  measured  on 
a  Cary  14  spectrophotometer.*  These  spectra  ranged  from  0.35  to  2.5  um,  an 
interval  practical  for  existing  remote  sensing  technology.  The  MOGS  algorithm 
was  evaluated  both  on  the  reflectance  spectra  from  which  it  was  derived  and  on 
spectra  collected  following  the  algorithm  development.  In  addition  digital  Images 
of  grain  size  distribution  and  moisture  distribution  were  developed  from  actual 
multispectral  scanner  data. 


2.  PROCEDURE 

The  research  involved  in  the  construction  of  the  MOGS  algorithm  was  divided 
into  two  parts.  First  the  Suits  radiative  transfer  vegetation  canopy  model  was 
converted  to  a  model  which  could  be  applied  to  beach  sands  (AQUASAND)  and,  second, 
the  MOGS  algorithm  was  developed  using  insight  obtained  from  the  AQUASAND  model. 

2.1  THE  AQUASAND  MODEL 

The  most  elementary  model  of  sand  reflectance  is  the  simple  plane  mixture 
model.  The  model  employs  the  assumptions  that  all  sand  particles  are  opaque 
and  are  randomly  mixed.  The  surface  of  the  sand  is  made  up  of  the  cross- 
sectional  areas  of  the  individual  particles.  The  model  calculates  the  bulk 
reflectance  as  being  a  weighted  average  of  the  particle  reflectances.  The 
model  fails  to  achieve  good  accuracy  because  the  transmittance  of  some  particles 
in  a  finely  divided  state  can  be  quite  large  and  multiple  scattering  between 
particles  is  often  significant. 

A  more  complex  model  and  the  one  used  in  this  work,  employs  concepts 
identical  to  those  employed  in  the  directional  reflectance  model  for  vegetative 
canopies  (Suits,  1972).  Although  a  vegetative  canopy  and  a  sand  profile  are 
visually  quite  different,  the  essence  of  the  reflection,  transmission,  and  multiple 
scattering  phenomena  is  much  the  same  for  both  cases. 

2.1.1  MODEL  ASSUMPTIONS 

The  first  assumption  of  the  model  is  that  the  scattering  components  are 
distributed  in  more  or  less  uniform  layers.  Because  of  the  wind  and  wave  action 
involved  in  the  creation  of  a  beach,  sand  profiles  tend  to  be  vertically  strati¬ 
fied  into  horizontal  layers  and  thus  the  assumption  is  satisfied. 

The  second  assumption  of  the  model  is  chat  spectral  flux  interacting  with 
the  sand  may  be  divided  into  two  types:  specular  and  diffuse.  Specular  flux 
represents  radiation  which  passas  between  Che  sand  grains  through  voids  and 
cracks  without  deviation.  All  incoming  radiation  is  specular  before  it  reaches 
the  sand  surface.  Diffuse  flux  is  some  fraction  of  the  specular  flux  which  has 
been  intercepted  by  a  scattering  component  (sand  grain)  at  least  once.  Following 
this  interception  it  may  be  scattered  forward  or  backwards ,  either  interacting 
with  ocher  sand  particles,  or  heading  upwards  to  the  sensor. 

A  third  assumption  of  the  model  is  that  Che  individual  particles  making  up 
the  sand  can  be  represented  by  a  sec  of  horizontal  and  vertical  Lambertian 
panels  which  have  the  same  reflectance  and  transmittance  properties  as  do  Che 
sand  particles.  This  assumption  defines  a  simplified  scattering  phase  function 
which  allows  the  calculation  of  ensemble  reflectances  in  closed  form.  Fortunately 


'(The  Cary  14  spectrophotometer  is  a  device  which  is  capable  of  digitally 
recording  the  reflectance  spectrum  of  a  surface  in  the  .35  to  2.5  .im  range.) 
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in  sands,  as  in  vegecacive  canopies,  multiple  scattering  dominates  over  phase 
function  scattering  because  of  the  high  density  of  the  scattering  components. 

This  characteristic  greatly  simplifies  the  necessary  model  calculations. 

The  fourth  assumption  of  che  model  is  that  the  diffuse  flux  moves  in  a 
Lambertian  manner,  both  upward  and  downward,  as  a  first  approximation.  The 
ensemble  reflectance  need  not  be  perfectly  Lambertian,  however,  it  should  be 
approximately  so.  Both  vegetation  and  sand  meet  this  criterion. 

As  mentioned  earlier  the  AQUASAND  model  operates  on  the  reflectance  (0) 
and  transmittance  (r)  of  a  sec  component  minerals  for  a  given  grain  size  to 
predict  che  reflectance  spectrum  of  a  beach.  One  method  of  obtaining  the 
necessary  o  and  t  values  would  be  to  measure  them  exactly  for  each  potential 
grain  size.  This  would  potentially  require  many  measurements  and  considerable 
sample  preparation.  Another  method,  and  the  one  used  in  this  study  would  be  to 
derive  the  components  of  transmittance  and  reflectance  for  each  mineral  and 
analytically  predict  che  o  and  t  for  any  given  grain  size.  The  necessary 
components  needed  to  derive  such  information  are  absorption  (a) ,  internal 
scattering  (s),  and  Che  forward  scattering  fraction  (FS) .  These  basic  measure¬ 
ments  were  obtained  by  measuring  che  t  and  0  of  3  chin  sections  of  each  potential 
mineral  type.  Each  thin  section  was  a  slightly  different  thickness  chan  che 
others.  By  using  an  iterative  curve  fitting  procedure  the  a,  s,  and  FS  values 
were  derived  using  che  Duntley  equations  (Duncley,  1942).  With  these  values 
calculated  at  10  nm  interval  throughout  the  spectrum  we  were  able  to  predict  the 
transmittance  and  reflectance  spectra  of  a  given  mineral  and  any  grain  size  we 
chose.  The  o  and  t  values  were  then  entered  into  the  AQUASAND  model  to  produce 
a  bulk  sand  reflectance  spectrum  (Figure  1). 

Moisture  content  was  entered  into  the  model  by  adding  an  appropriate  water 
spectral  absorption  factor  at  each  wavelength  for  which  a  spectral  value  was 
computed  for  “he  sand.  3y  entering  che  desired  percent-by-volume  water  concent 
the  correct  attenuation  due  to  water  was  computed  within  the  optical  pachway 
equations. 

2.1.2.  MODEL  RESULTS  AND  EVALUATION 

By  manipulating  the  mineralogical ,  grain  size,  and  moisture  model  input 
parameters  we  were  able  to  "create"  any  sand  type  we  wished.  Model  generated 
spectra  were  compared  to  empirical  spectra  obtained  from  seven  diverse  beach 
types  (Figure  2) .  The  empirical  measurements  were  made  on  a  Cary  14  spectro¬ 
photometer  which  allows  both  continuous  scanning  of  spectra  in  che  .35  to  2.5  Lm 
range  and  digital  recording  of  the  output  spectra. 

The  carbonate  beach  type  (Figure  2a)  is  composed  of  997.  exoskelatal  fragments 
of  marine  organisms  with  approximately  17.  dark  organic  debris  mixed  in.  Carbonate 
sands  are  characteristically  high  in  reflectance  with  high  reflectance  in  the 
"red"  spectral  region  (0.6  urn).  As  is  apparent  in  Figure  2  the  model  captures 
this  characteristic  quite  well. 

Figure  2b  and  c  are  spectra  of  predominately  iron-stained  quartz  beaches 
characteristic  of  the  Delaware  coastline.  In  general  these  spectra  are  depressed 
due  to  che  iron  scainlng  and  have  substantial  amounts  of  feldspar  incorporated 
with  chem. 

Figure  2d  is  a  non-iron  stained,  987.  pure  quartz  beach.  This  beach  type 
is  typical  of  chat  found  on  che  Gulf  of  Mexico  coastline.  These  beaches  have 
high  reflectance  due  to  che  lack  of  iron  stain  and  ocher  dark  minerals.  Besides 
quartz  there  are  trace  amounts  of  carbonate  and  organic  matter  amounting  to 
approximately  27.. 

Figure  2e  and  f  are  spectra  of  iron- stained  beaches  from  che  Lake  Michigan 
coastline.  Both  sands  exhibit  relatively  low  reflectance  although  spectrum  2f 
is  definitely  che  lower  of  the  two.  This  difference  is  related  to  the  more 
intense  iron-staining  and  larger  percentage  of  opaque  rock  fragments  in  the  latter 
beach  type. 
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Figure  2g  ip  a  heavy  mineral  beach  spectra  collected  on  the  coast  of  Oregon, 
which  exhibits  a  characteristically  low  reflectance  in  the  visible  region.  This 
is  due  to  a  high  amount  of  iron-staining  in  addition  to  a  large  percentage  of 
high  density,  opaque  minerals  (i.e.,  llmenite,  magnetite,  etc.).  Beaches  of 
this  type  tend  to  have  relatively  low  percentages  of  quartz,  on  the  order  of 
30  to  40  percent,  coupled  with  equal  amounts  of  feldspar.  Except  for  minor  errors, 
the  AQUASAND  model  correctly  predicted  the  spectrum  of  the  empirically  measured 
sand  in  all  3even  cases.  With  this  successful  validation,  work  began  on  the 
algorithm  predicting  mineralogy,  moisture,  and  grain  size  based  on  sand 
reflectance  spectra  in  the  .33  to  2.S  ym  range. 

2.2  DEVELOPMENT  OF  THE  MOGS  ALGORITHM 

From  Inspection  of  both  the  AQUASAND  generated,  and  empirical  spectra  it 
became  apparent  that  mineralogy  had.  by  far  the  greatest  influence  on  bulk 
reflectance.  So  great  is  this  influence  that  it  tends  to  mask  the  more  subtle 
features  of  changes  in  grain  size  and,  to  a  much  lesser  extent,  moisture.  It  . 
was  decided  that  in  order  for  the  algorithm  to  handle  a  broad  variety  of  sand 
types  and  still  maintain  the  resolution  needed  to  detect  small  spectral  features, 
a  preprocessing  clar sif ication  of  mineralogy  would  be  necessary. 

2.2.1  THE  CLASSIFICATION  OF  MINERALOGY 

In  order  to  discriminate  mineralogy  a  vector  length  decision  framework 
was  used.  The  concept  is  developed  as  follows. 

Suppose  that  there  are  two  points.  A  and  B,  located  in  two  dimensional 
space.  The  distance,  or  vector  length,  L,  from  A  and  B  can  be  expressed  in 
terms  of  the  X  and  Y  coordinate  locations  of  points  A  and  B  as: 

L  -  /  <XA  -  XB)2  ♦  (Ya  -  Yb)2  (1) 


This  is,  of  course,  related  to  the  Pythagorean  Theorem.  Now  suppose  we  have  a 
p-dimensional  system  with  A  and  B  located  in  each  dimension.  The  vector  length 
can  be  expressed  as 


L  - 


•  X1«>2. 


(2) 


where  is  che  location  of  point  A  in  the  ith  dimension  and  X.B  is  the 
location* of  point  B  in  the  ith  dimension. 


This  rationale  can  be  used  to  classify  some  point,  T,  as  being  the  member 
of  one  of  n  classes  (Aj ,  J»l,n),  by  finding  the  minimum  vector  length  from 
T  to  Aj(j-l.n).  In  other  words  T  is  said  to  be  a  member  of  the  class  which  is 
closest  to  it,  on  the  average,  across  all  p  dimensions.  The  minimum  vector 
length  is  defined  as 


min 


min 

. nl 


(X 


IT 


■  w 


(3) 


Notice  that  equation  3  has  no  provision  for  variability  in  the  n  classes,  there¬ 
fore,  L.  is  chosen  as  being  the  shortest  linear  vector  length.  If  each  class 
has  the  same  variability  associated  with  it  this  causes  no  difficulty.  In  this 
experiment,  however,  there  were  considerable  differences  in  variability  between 
the  classes  so  chat  a  modification  of  equation  3  had  to  be  made.  '  The  standard 
deviation  (SD)  was  used  to  modify  the  distance  between  T  and  Aj  related  to  each 
dimension  thus  removing  the  effects  of  variability  from  each  class.  This 
normalized  minimum  distance  equation  is  expressed  as  _ 


min 


■  min 

j-1 . n 
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(X 


it 
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(4) 


773 


where  SD,,  is  che  standard  deviation  associated  with  the  jctl  class  in  the  ith 
dimension. J 

In  che  application  of  this  method  to  the  classification  of  mineralogy,  the 
"dimensions"  are  spectral  bands  or  ratios  of  spectral  bands  and  che  "classes" 
are  mineralogical  types.  Eight  spectral  bands  (Table  1)  and  all  possible  unique 
ratios  of  chose  spectral  bands  were  used  to  classify  the  mineralogical  type  of 
an  input  sand  as  one  of  five  categories  (Table  2).  The  object  was  to  make  each 
category  as  homogeneous  as  possible  so  that  the  moisture  and  grain  size  regres¬ 
sions  which  followed  would  be  sensitive  to  small  scale  spectral  changes. 

2.2.2.  THE  DEVELOPMENT  OF  MOISTURE  AND  GRAIN  SIZE  REGRESSIONS 

Using  che  AQUASAND  model  we  found  that  information  related  to  moisture 
content  of  sands  is  best  derived  from  the  spectral  region  beyond  1.0  urn.  This 
is  due  to  che  fact  chat  che  spectral  reflectance  of  sand  in  this  region  is 
reduced  by  absorption  in  proportion  to  the  amount  of  water  present.  Exceptionally 
high  spectral  absorption  is  noced  near  1.4  and  1.9  urn.  Although  Che  spectral 
reflectance  in  these  regions  is  highly  correlated  to  moisture  we  did  not  consider 
them  since  acmsopheric  absorption  prohibits  their  use  by  an  airborne  sensor. 

Changes  in  grain  size  seem  to  manifest  themselves  most  clearly  in  che 
shorter  wavelengths  (.4-. 7  urn).  Grain  size  information  is  gained  by  light 
being  reflected  from  sand  grains  below  che  surface  through  surface  grains.  The 
'transmittance  through  che  surface  grains  is  reduced  by  internal  scattering  and 
absorption  of  the  particle.  Since  both  of  these  factors  are  dependent  on  thick¬ 
ness,  che  bulk  reflectance  of  a  sand  is  dependent  to  some  degree  on  che  grain 
size.  Theoretically  a  large  grain  sand  should  have  a  lower  reflectance  than  a 
small  grain  sand  of  che  same  mineralogical  composition  and  with  similar  moisture 
content.  According  to  our  measurements  this  appears  to  be  Che  case. 

This  grain  size  phenomenon  can  be  confounded  in  two  ways.  First  if  there 
is  no  scattering  or  absorption  within  che  grains  (i.e.,  a  perfectly  clear 
material  at  all  wavelengths)  there  can  be  no  attenuation.  Fortunately,  even  in 
our  purest  quartz  sands  there  were  enough  impurities  and  inclusions  to  give  some 
attenuation.  Second,  Che  sand  grains  may  be  opaque  and  thus  attenuate  coo  much 
light.  This  appears  to  be  the  case  in  the  heavy  mineral  and  carbonate  beaches. 

Most  of  the  bulk  reflectance  for  these  two  types  was  due  to  surface  reflectance  and 
little  if  any  was  due  to  light  transmitted  through  che  surface  grains  from  below. 

We  were  unable  to  create  accurate  grain  size  equations  for  these  types. 

Utilizing  che  physical  phenomena  discussed  above  we  were  able  to 
develop  multiple  linear  regression  equations  for  predicting  moisture  in  all 
five  mineralogical  classes  and  grain  size  for  three  of  che  five  mineralogical 
classes.  The  basis  for  all  che  regressions  except  one  was  che  sample  group 
corresponding  to  a  given  mineralogical  class.  The  single  exception  was  che 
grain  size  equation  corresponding  to  a  pure  quartz  beach.  Our  sampi_(-  within 
this  type  consisted  of  a  single  grain  size  (0.22  mm)  and,  as  such  did  not  provide 
an  adequate  basis  for  regression  equations.  For  this  case  we  used  AQUASAND 
generated  spectra  to  simulate  a  wide  range  of  grain  sizes  in  order  to  add  grain 
size  variability  to  che  data  set. 

Seventeen  spectral  bands  between  0.4  and  2.5  urn  were  chosen  for  use  in  che 
regressions  (Table  3).  Within  che  17  bands,  only  chose  which  were  predicted  by 
che  AQUASAND  model  to  be  most  informative  were  used.  In  this  way  we  could  be 
reasonably  certain  chat  the  regression  equations  would  respond  the  correct 
parameter  and  thus  yield  accurate  predictions.  The  predictive  equations 
together  with  che  associated  standard  errors  (SE),  and  coefficient  of  variation 
(RZ)  are  given 'in  Table  4. 

In  summary  the  MOGS  algorithm  {Figure  3)  represents  a  compucer  controlled 
package  of  equations.  The  input  is  a  sec  of  17  spectral  reflectance  bands 
obtained  from  an  unknown  sand.  Based  on  these  bands,  the  sand  is  classified  as 
being  a  member  one  of  five  mineralogical  types.  Depending  on  che  mineralogical 
type,  che  appropriate  moisture  and  grain  size  (where  applicable)  equations  are 
applied  to  the  data.  The  output  from  che  MOGS  algorithm  is  che  predicted  mineralogi¬ 
cal  class,  che  predicted  moisture,  and  che  predicted  grain  size. 
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3.  TEST  RESULTS  OF  THE  MOGS  ALGORITHM  OH  LABORATORY  SPECTRA 


The  MOGS  algorithm  was  first  tested  on  70  of  the  81  samples  from  which  it  was 
derived  and  che  results  were  very  promising  as  can  be  seen  in  Table  5.  The 
classification  of  mineralogy  was  99+7.  correct.  The  overall  correlation  of 
predicted  to  actual  moisture  was  96 X  (significant  at  the  .001  level)  and  the 
overall  correlation  of  predicted  to  actual  grain  size  was  887.  (significant  at 
the  .001  level).  However,  testing  any  equation  or  algorithm  on  the  samples 
from  which  it  was  derived  is  not  conclusive.  For  this  reason  the  MOGS  algorithm 
was  tested  on  several  ocher  beach  sand  samples  which  were  independently  collected 
and  spectrally  measured  following  the  algorithm  construction.  These  results 
ire  given  in  Table  6.  In  each  case  the  MOGS  algorithm  selected  a  mineralogy 
which  allowed  the  moisture  and  grain  size  regression  to  operace  correccly 
The  independent  test  yielded  an  actual  moisture  to  predicted  moisture  correlation 
of  .95  (significant  at  the  .01  level).  The  prediction  of  grain  size  was  in  no 
case  more  than  0.07  mm  different  from  the  actual  grain  size. 

4.  .TEST  OF  THE  MOGS  ALGORITHMS  ON  MSS  DATA 

The  next  logical  test  for  che  MOGS  algorithm  was  to  evaluate  it  on  actual 
MSS  data.  Such  an  investigation  is  presently  underway  using  data  obtained  from 
cne  Environmental  Research  Institute  of  Michigan  (ERIM)  multispectral- scanner 
(MSS)  flown  over  a  portion  of  the  Lake  Michigan  shoreline  (November  1,  1978). 

The  MOGS  algorithm  was  modified  to  conform  to  che  12  band  configuration  of  che 
ERIM  scanner. 


Although  a  complete  discussion  or  Che  Lake  Michigan  MSS  test  is  beyond  che 
scope  of  this  paper,  preliminary  results  look  quite  good  (Table  7).  The  correla¬ 
tion  of  predicted  to  actual  moisture  concent  is  .91  (significant  at  che  .01  level) 
and  the  prediction  of  grain  size  is,  in  no  case,  greater  chan  .09  mm  different 
from  che  actual  grain  size. 

The  moisture  prediction  appears  particularly  poor  for  high  moisture  concents. 
This  may  be  due  to  che  fact  chat  wet  sand  at  the  test  sites  appeared  to  exhibit 
bi-directional  dependencies  (t.e. ,  a  failure  to  behave  in  a  Lambertian  manner). 
These  bi-directional  characteristics  are  enhanced  at  low  sun  angles  and  are  not 
accounted  for  by  the  MOGS  algorithm.  Although  che  flight  cook  place  at  1:30  EST 
the  sun  was  only  38  above  the  horizon  on  November  1.  To  minimize  bi-directional 
reflectance,  future  aircraft  flights  should  be  made  during  complete  mid-altitude 
(3000  m)  cloud  cover  or  sunny  skies  with  Che  sun  close  to  che  zenith  (summer  sun) . 

Applying  che  same  moisture  and  grain  size  equations  used  in  the  previous 
analysis,  che  entire  Pentwacer  State  Park  beach  on  a  pixel  by  pixel  basis  (in 
this  case  1.5  x  1.5  meters)  was  classified  with  respect  to  grain  size  and  moisture 
concent.  The  two  MOGS  generated  digital  maps  (see  Figure  4)  show  the  predicted 
moisture  and  grain  size  distributions  on  che  beach  at  Pentwacer.  The  ground  truth 
measurements  taken  at  the  clme  of  flight  correlate  well  with  these  images. 

Figure  4  helps  to  demonstrate  how  an  encire  sandy  coastline  could  be  analyzed  in 
respect  to  moisture,  grain  size,  and  gross  mineralogy  using  a  smell  subsection 
as  calibration. 


5.  CONCLUSION 

The  development  of  che  MOGS  algorithm  has  demonstrated  che  feasibility  of 
obtaining  quantitative  moisture  and  grain  size  information  from  che  spectral 
reflectance  of  beach  sands.  The  determination  of  grain  size  is  dependent  on  the 
sand  grains  being  neither  opaque  or  perfectly  clear. 

The  two  stage  nature  of  che  MOGS  algorithm  is  directly  responsible  for  its 
broad  applicability  without  loss  of  detail.  By  separating  che  mineraloglcal 
types  prior  to  che  prediction  of  moisture  and  grain  size  much  of  che  variability 
which  could  easily  hide  small  scale  changes  is  removed.  The  use  of  a  vector 
length  discriminant  function  to  classify  mineralogy  worked  extremely 
well  in  this  application,  since  36  different  dimensions  could  be  simultaneously 
evaluated.  The  use  of  a  multistage  approach  involving  multiple  classification 
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cecnniques  is  a  powerful  cool;  one  which  will  very  likely  be  useful  in  many 
areas  of  remoce  sensing. 

The  Lake  Michigan  field  test  has  furcher  demonstrated  che  MOGS  algorithm's 
applicability  to  actual  remotely  sensed  data.  Grain  size  was  predicted  to  within  09 
mean  diameter  of  actual  while  beach  moistures  less  than  20%  were  accurately  pre- 
dieted.  In  all  cases  the  computer  algorithm  correctly  identified  che  Michigan 
beach  mineralogy  as  being  a  predominantly  quartz  iron  stained- feldspar  beach. 


Table  1.  The  S  spectral  bands  used  in  che  breakdown  of 
beach  mineralogy  into  1  of  5  categories.  In 
addition  to  these  8  bands  all  unique  ratio 
combinations  were  also  used. 


Band  4 

2 

3 

4 

5 

6 

7 

8 


Wavelength  Range 

- .  43-  47 - 

.47-. 49 
.51-. 53 
.53-. 56 
.59-. 63 
.30-  .90 
.90-1.0 
1.0-1. 1 


LdSl 


Table  2.  The  5  potential  mineralogical  classifications 


Class 

1  - 

2 

3 

4 

5 


Description 

Iron  stained  Atlantic  coast  type 
Iron  stained  Michigan  coast  type 
Iron  stained  pure  quartz  type 
Heavy  mineral  type 
Carbonate  type 


Table  3.  The  17  bands  used  in  the  development  of  moisture 
and  grain  size  regression  equations 


Band  # 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


Wavelength  Range  (urn) 

- 0. 53-OV47 

0.47-0.49 

0.49-0.51 

0.51-0.53 

0.53-0.56 

0.56-0.59 

0.59-0.63 

0.63-0.67 

0.70-0.75 

0.75-0.80 

0.80-0.90 

0.90-1.00 

1.00-1.10 

1.10-1.20 

1.20-1.35 

1.50-1.85 

2.10-2.50 
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Table  4.  Multiple  linear  regression  equations  for  the 
prediction  of  moisture  and  grain  size.  The 
equations  are  listed  by  mineralogical  class. 

Grain  size  is  in  mm. 

a.  Iron  stained  quartz  -  Atlantic  coast 

Predicted  moisture  %  -  67.964  -  65.046  [%) 

S.E.  -  3.08%.  R2  -  0.388 

Predicted  grain  size  -  6.87  -  3.4634  (Band  7)I/4  +  .0300  (Band  1)  +  ,01672(Band  15) 
S.E.  -  0.13  mm.  R2  -  0.603 

b.  Iron  stained  quartz  -  Michigan  coast 

Predicted  moisture  %  -  60.149-49.961  (|gj|  ffi-2.226  [7) 

S.E.  -  2.56%.  R2  -  .970 

Predicted  grain  size  •  0.6405-0.0152  (Band  5)  -  .0047  (Band  17) 

S.E.  -  0.055  m.  R2  -  0.558 

c.  don- Iron  stained  quartz 

Predicted  Moisture  %  -  127.02-65.159  (jj~j~^§)  -  64.054  (||2£[_£|) 

S.E.  -  2.12%,  R2  -  0.971 

Predicted  grain  size  -  1.158-2.328  (Band  10)  +  .3201  (g-jjg|  [)  +  0.2858  (Band  Id 
S  E.  and  R2  not  applicable 

d.  Carbonate 

Predicted  moisture  %  -  596.28-642  (|^jj- T7>~  1-081  (Band  14>  +  °-^38  (Band  17) 

S.E.  -  4.09%.  R2  -  .879  No  grain  size  equation. 

e .  Heavy  mineral 

Predicted  Moisture  %  -  19.284+11.194  (§gB|- |j)  -  1.081  (Band  14)  +  0.1538  (Band  17) 
S.E.  *  4.09%,  R2  ■  .879  No  grain  size  equation. 
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Table  5. 


Sand  I,D. 

Actual  Predicted 


A1  A 

A2  A 

A3  A 

A4  A 

A5  A 

A6  A 

A7  A 

Afi  A 

A9  A 

A10  A 

31  B 

B2  3 

33  B 

B4  B 

33  B 

B6  B 

87  B 

B6  B 

B9  B 

BIO  B 

Bll  B 

B12  B 

B13  A 

314  B 

B15  B 

B16  B 

B17  B 

B18  B 

E19  B 

B20  B 

Ml  M 

M2  M 

M3  M 

M4  M 

M5  M 

M6  M 

M7  M 

M8  M 

M9  M 

MIO  M 

Mil  M 

M12  M 

M13  M 

M14  M 

M15  M 

M16  M 

01  H 

02  H 

03  H 

04  H 

05  H 

06  H 

07  H 


Comparison  of  actual  sand  parameters  to  predicted 
classification  by  the  MOGS  algorithm 


Moisture  7. 


Actual 

Predicted 

4.5 

3.8 

29.4 

26.4 

15.0 

11.8 

28.4 

28.0 

11.3 

12.9 

33.4 

27.1 

13.1 

13.8 

24.9 

27.3 

8.8 

19.4 

29.7 

27.9 

21.0 

17.3 

24.6 

25.5 

14.2 

15.1 

27.2 

26.3 

'11.0 

8.9 

19.0 

25.6 

26.8 

28.2 

24.9 

21.0 

6.2 

5.9 

★ 

7.5 

20.0 

23.0 

34.0 

33.1 

6.0 

4.4 

31.0 

26.0 

18.0 

22.2 

32.0 

31.5 

18.0 

20.8 

23.0 

23.9 

3.0 

7.1 

23.0 

27.8 

5 

8.2 

15 

13.4 

25 

26.2 

30 

27.4 

3 

5.5 

15 

15  7 

0 

0.0 

0 

0.0 

0 

0.0 

0 

0.0 

10 

12.0 

25 

24.4 

30 

27.2 

10 

12.8 

20 

20.6 

15 

16.8 

0 

0 

5 

6.4 

10 

11.7 

15 

20.3 

20 

22.0 

25 

21.9 

30 

27.5 

Grain  Size  mn 


Accual 

Predicted 

.35 

.36 

.35 

.46 

.37 

.40 

.50 

.48 

.43 

49 

.32 

.39 

.35 

.43 

.38 

37 

.  44 

.30 

.43 

.44 

.40 

.62 

.76 

.58 

.46 

.57 

.88 

.76 

.63 

.32 

95 

.76 

.  71 

.76 

.76 

.63 

.71 

.62 

.71 

.57 

.81 

.76 

.69 

.73 

.55 

.60 

.83 

.71 

.67 

.69 

.57 

.72 

.56 

.67 

.65 

.76 

.94 

.68 

.60 

.61 

.26 

.27 

.36 

.20 

.36 

.34 

.36 

.36 

.41 

.38 

.41 

.45 

.23 

.31 

.29 

.29 

.41 

.39 

.36 

.35 

.23 

.31 

.13 

.28 

.23 

.25 

.41 

.36 

.41 

.40 

.23 

.26 

no 

grain 

size 

no 

grain 

size 

no 

grain 

size 

no 

grain 

size 

no 

grain 

size 

no 

grain 

size 

no 

grain 

size 

*No  moisture  data 
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Sand 

I.D. 

Table  5.  (Continued) 
Moisture  7. 

Grain  Size  mn 

Actual 

Predicted 

Actual 

Predicted 

Actual 

Predicted 

08 

H 

35 

24.7 

no 

grain 

size 

MX1 

MX 

25 

24.6 

.22 

.  32 

MX2 

MX 

0 

0.0 

.22 

.  10 

MX3 

MX 

10 

14.1 

.22 

.  13 

MX4 

MX 

15 

18.1 

.22 

.22 

MX5 

MX 

30 

34.2 

.22 

4 

MX6 

MX 

35 

33.7 

.22 

34 

MX7 

MX 

20 

23.1 

.22 

.33 

MX8 

MX 

5 

6.9 

.22 

.11 

Cl 

C 

0 

0 

no 

grain 

size 

C2 

c 

10 

7.8 

no 

grain 

size 

C3 

c 

5 

5.7 

no 

grain 

size 

C4 

c 

20 

13.4 

no 

grain 

size 

C5  . 

c 

15 

18.8 

no 

grain 

size 

C6 

c 

30 

30.8 

no 

grain 

size 

C7 

c 

25 

27.9 

no 

grain 

size 

C8 

c 

40 

38.4 

no 

grain 

size 

C9 

c 

50 

46 . 4 

no 

grain 

size 

Table  6.  Comparison  of  actual  parameters  to  predicted 

classif icattions  by  tne  MOGS  algorithm.  Samples 
used  here  were  collected  independently  of  chose 
on  which  Che  algorithm  is  based. 


Sample 

Mineralogy 

Moisture 

7. 

Grain  Size 

(mm) 

Actual 

Predicted 

Actual 

Predicted 

Actual 

Predicted 

MlCHl" 

"  MICH 

37.7' 

25. 6 

.26  "■ 

.23 

MICH2 

MICH 

3.2 

3.7 

.28 

.30 

MICH3 

MICH 

0.3 

0.0 

.25 

.32 

MICH4 

MICH 

12.1 

16.0 

.23 

.29 

MICH5 

MICH 

15.0 

12.7 

.25 

.30 

MICH6 

MICH 

28.0 

26.2 

.28 

.29 

Table  7.  Comparison  of  actual  parameters  co  predicted 

classifications  by  the  MOGS  algorithm.  The  sand 
spectra  used  here  were  collected  by  the  ERIM  MSS. 


Sample 

Mineralogy 

Moisture 

7. 

Grain  Size 

(mm) 

Actual 

Predicted 

Actual 

Tredicted 

Actual 

Predicted 

MICH1 

MICH 

22.3 

12. 1 

.25 

.32 

MICH  2 

MICH 

1.0 

0.0 

.23 

.28 

MICH3 

MICH 

8.0 

10.1 

.22 

.25 

MICH4 

MICH 

28.0 

16.9 

.26 

.35 

MICH5 

MICH 

5.0 

2.2 

.24 

.26 
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Figure  L.  Flow  diagram  of  the  AQUASAND  MODEL. 
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OUTPUT  INFORMATION 


FIGURE  3.  FLOW  DIAGRAM  OF  THE  MOGS  ALGORITHM. 
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MOISTURE  DISTRIBUTION  IMAGE  -  LIGHTER  SHADE  INDICATES  HIGHER  MOISTURE 

CONTENT 


GRAIN  SIZE  DISTRIBUTION  LIAGl  -  LIGHTER  SHADE  INDICATES  LARCER  GRAIN  SIZE 
.  15-  .  25  mm  Hi  .  25  -  .  AO  mm  Hi  ■  AO  mm 


PANCHROMATIC  AERIAL  PHOTOGRAPH 

FIGURE  A.  DIGITAL  IMAGERY,  GENERATED  BY  THE  MOGS  ALGORITHM  SHOWING  THE  DISTRIBUTION 
OF  MOISTURE  AND  GRAIN  SIZE  ON  PENTWATER  BEACH  (PENTWATER  STATE  PARK 
MIChlGAu) .  WHITE  AREAS  ARE  EITHER  UNCLASSIFIED  REGIONS  OR  OPEN  W/ "ER 
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ABSTRACT 


An  algorithm  was  developed  which  determines  the 
mineralogy,  moisture,  and  grain  size  of  beach  sands 
based  on  the  hemispherical  reflectance  in  17  discrete 
spectral  bands.  The  bands  chosen  range  between  0.40 
and  2.5  pm,  a  wavelength  range  practical  for  existing 
multispectral  remote  sensing  technology.  The  sand 
rpectra  on  which  the  mineralogy,  moisture,  and  grain 
size  algorithm  (M0GS)  is  based  were  obtained  from 
laboratory  spectrophotometric  measurements.  Selected 
spectral  bands  are  used  in  a  vector-length-decision 
framework  to  determine  the  mineralogical  class  of  the 
input  sand.  Multiple  linear  regressions  are  then 
used,  within  a  given  mineralogical  class,  to  determine 
the  moisture  and  grain  size  of  the  sand. 

The  predictive  results  of  the  M0GS  algorithm  are 
very  encouraging.  When  tested  on  70  of  the  sand 
reflectance  spectra  from  which  it  was  derived,  the 
correlation  of  actual  to  predicted  moisture  and  grain 
size  was  96%  and  88%,  respectively. 

The  M0GS  algorithm  has  been  successfully  tested 
using  aircraft  multispectral  scanner  (MSS)  data 
collected  over  the  Lake  Michigan  shoreline.  The 
algorithm  correctly  identified  gross  mineralogy  and 
predicted  grain  size  to  within  0.09  mm  of  nteasui  od 
values.  Some  difficulties  were  encountered  in 
predicting  high  beach  sand  moistures,  probably  due 
to  the  increasing  non-Lambertian  nature  of  sand  as 
the  moisture  concent  of  the  sand  increased. 


INTRODUCTION 


During  the  past  ten  years  remote  sensing  has  been  proven  capable  of 
delineating  outcrop  lithologies,  a  procedure  often  carried  out  in  the  early 
stages  of  mineral  exploration.  Typically,  remote-sensing  spectral  data  are 
compared  to  those  in  a  reference  library  containing  laboratory-acquired 
spectral  data  of  individual  minerals  and  rocks.  Through  the  use  of  a 
computer,  correlations  are  made  between  the  reference  and  remote-sensing 
spectra  in  an  attempt  to  identify  the  surface  material  (Vincent,  Thomson 
and  Watson,  1972) . 

In  this,  investigation, -beach  sands  were  analyzed  with  the  intent  of 
determining  not  only  mineralogy  but  also  moisture  and  grain  size.  These 
three  parameters  are  of  interest  from  both  „  "ological  and  engineering  points 
of  view.  The  mean  grain  size  at  a  reference  point  on  a  beach  is  a  funda¬ 
mental  characteristic  of  the  beach.  Studies  have  indicated  that  mean  grain 
size  is  related  to  beach  face  slope  (Komar,  1976),  water  percolation  and 
permeability  and  ease  of  grain  movement  (Komar,  1977;  Zenkovich,  1967; 

Huntley  and  Bowen,  1975;  Fraser  and  Hester;  1977;  Self,  1977).  These  factors 
play  a  majc  .  role  In  determining  if  a  beach  will  be  erosional  or  depositional 
under  different  wave  conditions  (Madsen  and  Grant,  1976;  Swart,  1976). 

Additionally,  recognition  of  these  three  sand  parameters  would  also 
allow  the  identification  of  beach  mineral  deposits  based  on  grain  size  and 
mineralogy.  Since  beaches  are  formed  by  intense  erosion  of  the  parent 
materials,  more  resistant  minerals  tend  to  be  preserved  and  concentrated 
while  others  are  dissipated.  Depending  on  the  origin  of  the  sand,  these 
residual  minerals  may  be  of  economic  value. 


1 


In  order  to  determine  the  mineralogy,  moisture,  and  grain  size  of  a 
given  sand,  a  two-stage  vector-leingth-decision  and  multiple  linear 
regression  approach  was  used.  This  entailed,  first,  the  breakdown  of 
mineralogy  into  five  categories.  Each  category  was  then  analyzed  for 
moisture  content  and  grain  size  using  predictive  multiple  linear  regressions 
based  on  selected  spectral  bands. 

BACKGROUND 

The  theoretical  basis  for  the  work  presented  in  this  paper  was  developed 
over  a  period  of  approximately  ten  years.  Studies  have  been  made  towards 
understanding  the  effects  that  grain  size,  moisture  coating,  and  mineralogy 
of  particulates  have  on  reflected  radiation  from  beach  sands.  It  was  shown 
by  Emslie  (1966)  and  Aronson,  et  al,  (1967)  that  reflected  radiation  from 
beach  sand  is  a  function  of:  the  wavelength  of  the  radiation;  the  optical 
constants  of  the  medium,  i.e.,  n  (refractive  index)  and  k  (index  of 
absorption);  the  particulate  grain  size;  the  packing  density;  and  the  rough¬ 
ness  of  the  surface.  Vincent  and  Hunt  (1968)  used  reflectance  data  based 
on  ground  laboratory  samples  to  improve  on  the  work  of  Aronson,  et  al.  The 
model  of  Vincent  and  Hunt  also  accounted  for  volume  scattering  of  reflected 
radiation  in  sand  particles. 

Leu  (1977)  working  at  the  Environmental  Research  Institute  of  Michigan 
(ERTM)  collected  beach  sands  (from  Delaware)  and  measured  their  reflectance 
ia  the  laboratory  in  the  0.4  to  2.5  p;i  range  using  a  Cary  14  laboratory 
spectrophotometer.  (The  Car-/  14  spectrophotometer  is  a  device  which  is 
capable  of  digitally  recording  the  reflectance  spectrum  of  a  surface  in  the 
0.35  to  2.5  um  range.)  Leu  then  used  multiple  linear  regression  to  correlate 


reflectance  values  In  discrete  reflectance  bands  to  moisture  and  grain  size. 
Leu's  work  showed  promise,  but  his  regressions  did  not  work  when  beach  sand 
reflectance  values  ocher  than  those  from  Delaware  beaches  were  used.  This 
indicated  the  mineralogy  of  a  beach  sand  strongly  affects  its  reflectance 
spectra  and  chat  the  effects  of  mineralogy  must  be  understood  before 
moisture  and  grain  size  prediction  algorithms  can  be  run. 

In  order  to  understand  fully  the  effects  of  changes  in  mineralogy, 
moisture,  and  grain  size  on  the  spectral  reflectance  of  beach  sand,  a  sand 
reflectance  model  was  used.  The  model,  an  adaptation  of  the  Suits  radiative 
transfer  vegetation  canopy  model,  is  known  as  AQUASAND  (Suits,  1972).  It 
uses,  as  inputs,  the  reflectance  and  transmittance  for  each  mineral 
comprising  the  beach  sand  (i.e.,  quartz,  feldspar,  magnetite,  etc.).  In 
addition,  input  to  the  model  includes  the  sand  grain  size,  the  void  space 
in  the  sand,  and  the  moisture  profile  as  a  function  of  depth.  By  varying 
these  input  parameters  insight  was  gained  into  the  effect  of  physical 
changes  on  the  bulk  sand  reflectance. 

Using  the  information  obtained  from  AQUASAND  (Shuchman,  et  al.,  1978) 
the  mineralogy,  moisture,  and  grain  size  (HOGS)  algorithm  was  developed 
using  reflectance  spectra  measured  on  a  Cary  14  spectrophotometer.  These 
spectra  ranged  from  0.35  to  2.5  yra,  an  interval  practical  for  existing  remote 
sensing  technology.  The  MOGS  algorithm  was  evaluated  both  on  the  reflectance 
spectra  from  which  it  was  derived  and  on  spectra  collected  following  the 
algorithm  development.  In  addition,  digita1  images  of  grain  size  distri¬ 
bution  were  developed  from  actual  multispectral  scanner  data. 


DEVELOPMENT  OF  THE  MOGS  ALGORITHM 


The  procedure  of  developing  an  algorithm  to  predict  mineralogy,  moisture 
and  grain  size  of  beach  sands  was  divided  into  three  segments.  First,  a 
data  base  of  sand  reflectance  values  was  selected  from  which  to  build  the 
algorithm,  second,  the  necessary  equations  were  developed,  and  third,  the 
algorithm  was  evaluated  using  actual  remotely  sensed  aircraft  data. 

All  of  the  equations  that  make  up  the  mineralogy,  moisture,  and  grain 
size  (MOGS)  algorithm  are  based  ultimately  on  81  laboratory  measured 
reflectance  spectra  of  beach  samples  obtained  from  a  Cary  14  spectrophoto¬ 
meter  operated  by  the  Environmental  Research  Institute  of  Michigan  (ER1M) . 

The  sand  samples  were  collected  from  five  diverse  beach  types  located  in 
various  coastal  areas  of  the  continental  United  States  between  June  1974 
and  October  1978.  The  use  of  a  large  range  of  sand  types  was  deemed 
necessary  to  give  the  MOGS  algorithm  a  wide  field  of  applicability.  The 
mean  grain  si^a  (diameter),  moisture  content,  and  approximate  location  of 
each  beach  sand  sample  used  in  the  development  of  the  MOGS  algorithm  are 
given  in  Table  1. 

From  inspection  of  both  the  AQUASAND  generated,  and  empirical  spectra 
it  became  apparent  that  mineralogy  had  by  far  the  greatest  influence  on  the 
reflectance  spectra.  So  great  is  this  influence  that  it  tends  to  mask  the 
more  subtle  features  of  changes  in  grain  size  and,  to  a  much  lesser  extent, 
moisture. 

In  order  to  achieve  the  fine  detail  needed  for  the  classification  of 
moisture  and  grain  size  while  still  maintaining  the  applicability  to  a 
large  range  of  mineralogies  a  two-stage  procedure  was  established  in 
the  development  of  the  MOGS  algorithm  (Figure  1).  The  first  stage  entails 


TABLE  1 


The  mean  grain  sizes  and  moisture  contents  of  the  81  sand  samples 
used  in  the  development  of  the  MOGS  algorithm.  Samples  are  given 

by  mineralogical  class. 


SanJ  I.D. 

Moir-curc 
Can:«\-.t  (*) 

Al 

4.5 

A2 

29.4 

A3 

15.0 

A4 

23.4 

>3 

1 1.3 

A6 

33.4 

A7 

13.1 

A8 

24.9 

A9 

8.S 

AIO 

29.7 

ei 

21.0 

B2 

24.6 

B3 

14.2 

F.'. 

27.2 

&j 

11.3 

to 

19.0 

n 

6.3 

E6 

24.9 

f>9 

6.2 

no 

21.3 

FIX 

20.0 

412 

34.0 

£13 

6.0 

B14 

31.0 

BIS 

13.0 

B16 

32.0 

E17 

13.0 

E1S 

23.0 

Big 

3.0 

B20 

23.0 

m 

5.0 

M2 

15.0 

M3 

25.0 

k; 

30.0 

MS 

5.0 

M6 

15.0 

117 

0.0 

M3 

0.0 

M9 

0.0 

MIO 

0.0 

Mil 

10.0 

Mi'an  Cf.’iin 


SI  sc  (isO 

Sand  I.D 

.35 

M12 

.35 

H13 

.37 

Ml  4 

.50 

Ml  5 

.43 

Ml  6 

.32 

1117 

.35 

MIS 

.33 

M19 

.44 

HI 

.43 

H2 

.40 

K3 

.76 

K4 

.46 

H5 

.38 

116 

.63 

H7 

.95 

Hi 

.71 

H9 

.76 

K 10 

.71 

HI  1 

.71 

1112 

.31 

MSI 

.69 

MS  2 

.55 

MS  3 

.83 

MX4 

.67 

1X5 

.57 

MX  6 

.56 

1X7 

.65 

MS  f, 

.94 

MS  9 

.60 

IK  10 

.36 

IK  ll 

.36 

MX  12 

.36 

Cl 

.36 

C2 

.41 

C3 

.41 

C4 

.23 

C5 

.29 

Co 

.41 

C7 

.36 

C3 

.23 

C9 

Moisture 

Moan  Crain 

Content  (2) 

Sire  (r.n) 

25.0 

.23 

30.0 

.23 

10.0 

.41 

20.0 

.41 

15.0 

.23 

22.0 

.23 

2.0 

.31 

8.0 

.26 

0.0 

.32 

5.0 

.32 

10.0 

.32 

15.0 

.32 

20.0 

.32 

25.0 

.40 

30.0 

.40 

35.0 

.40 

0.0 

.17 

15.0 

.22 

10.0 

.32 

5.0 

.40 

25.0 

.22 

0.0 

.22 

10.0 

15.0 

.22 

30.0 

.22 

35.0 

.22 

20.0 

.22 

5.0 

.22 

16.0 

.22 

4.0 

.22 

10.0 

.22 

29.0 

.22 

0.0  .  spprox 

lr.lx.25  wa 

10.0 

lxlx.25 

5.0 

lxlx.25 

20.0 

lxlx.25 

15.0 

Ixlx. 23 

30.0 

lxlx.25 

25.0 

lxlx.25 

40.0 

lxlx.25 

50.0 

lxlx.25 

here: 

A  *  Indian  River  Inlet,  Delaware 
B  ■  Dei^are  Bay,  Delaware 

M  ■  Michigan  Coastline  (Sleeping  Bear  State  Park,  Petosky  State  Park, 
Mason-Oceana  County  Line,  Pentwater  State  Park,  and  Muskegon 
State  Park) 

H  ■  South  Beach  and  Glen  Eden  3each,  Oregon  Coastline 
MX  *  Panama  City,  Florida  (Gulf  of  Mexico) 

C  *  Marine  Carbonate,  Florida  Keys 
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INPUT 


FIGURE  1.  FLOW  DIAGRAM  OF  THE  MOGS  ALGORITHM. 


a  breakdown  of  sand  mineralogy  into  discrete  classes.  Within  each  class, 
second  stage  multiple  linear  regressions  were  used  to  derive  the  moisture 
and  grain  size  information.  Using  this  method  the  second  stage  regressions 
do  not  have  to  account  for  the  large  spectral  effects  of  mineralogy  charac¬ 
teristic  of  diverse  beach  sands.  This  division  increased  the  accuracy  and 
predictive  ability  of  the  moisture  and  grain  size  regression  equations 
considerably. 

In  order  to  determine  which  spectral  regions  would  best  be  able  to 
differentiate  and  predict  the  parameters  of  interest,  the  AQUASAND  beach 
sand  model  was  used.  By  varying  the  model  input  parameters  in  a  logical 
fashion  it  was  possible  to  predict  what  regions  of  the  spectrum  yielded  the 
most  useful  information.  These  spectral  regions  were  then  used  in  the  MOGS 
algorithm  to  predict  the  parameters  of  interest.  Because  the  HOGS  algorithm 
is  to  be  used  on  aircraft  collected  multispectral  scanner  (MSS)  data,  all 
the  input  spectral  data  were  divided  into  17  spectral  "bands"  which  are 
feasible  (dictated  by  atmospheric  transmittance)  for  implenentation  using 
existing  MSS  technology  (Table  2) .  In  any  given  portion  of  the  MDGS  algorithm 
only  a  subset  of  these  17  bands  were  used,  as  will  be  shown  later. 

The  Prediction  of  Mineralogy 

Rather  than  attempting  to  predict  the  individual  mineral  components  of 
the  different  sands  it  was  decided  that  grouping  the  sands  into  homogeneous 
types  would  be  more  productive.  As  such,  five  groups  were  defined  which, 
for  the  most  part,  corresponded  with  the  geographic  location  of  collection 
of  the  data  base  samples.  The  five  mineralogical  classes  are  as  follows: 

1.  Iron-stained  Atlantic  coast  type  (A,  B) , 

2.  Iron-stained  Michigan  type  (M) , 

3.  Non-Iron-stained  pure  quartz  type  (MX), 
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TABLE  2 


The  17  spectral  bands  used  in  the  development  of  ! 
moisture  and  grain  size  regression  equations 


Band  If 

2 

3 

4 

5 

6 
7 

a 

9 

10 

11 

12 

13 

14 

15 

16 
17 


Wavelcncch  Ranee  (um) 

- 07TOJ747 

0.47-0.49 
0.49-0.51 
0.51-0.53 
0.53-0.56 
0.56-0. 59 
0.59-0.63 
0.63-0.67 
0.70-0. 75 
0.75-0.30 
0.80-0.90 
0.90-1.00 
1.00-1.10 
1.10-1.20 
1.20-1.35 
1.50-1.85 
2.10-2.50 


TABLE  3 

The  eight  spectral  bands  used  in  the  breakdown  of 
beach  mineralogy  into  one  of  five  categories.  In 
addition  to  these  eight  bands  all  unique  ratio 
combinations  were  also  used. 


Band  If 

— r-^ 
2 

3 

4 

5 

6 

7 

8 


Wave  length  Ranee  (ur.i) 

- TZttT - 

.47-. 49 
.51-. 53 
.53-. 56 
.59-. 63 
.80-. 90 
.90-1.0 
1.0-1. 1 


4.  Heavy  mineral  (dark  sand)  type  (H) 

5.  Carbonate  type  (C). 


Typical  spectra  of  these  five  miner alogical  categories  is  given  in 
Shuchman,  et  al,  (1978).  To  discriminate  the  five  mineralogical  classes 
a  minimun- vector-  length  decision  framework  was  used.  The  concept  is 
developed  as  follows.  Suppose  that  there  are  two  points,  A  and  B,  located 
in  two-dimensional  space.  The  distance,  or  vector  length,  L,  from  A  to  B 
can  be  expressed  in  terms  of  the  X  and  Y  locations  of  points  A  and  B  as 

l  -  [<xA  -  x/  +  <ya  -  y2]  %  .  (1) 

This  is,  of  course,  related  to  the  Pythagorean  Theorem.  Now  suppose  we 
have  a  p-dimensional  system  with  A  and  B  located  in  each  dimension.  The 
vector  length  can  be  expressed  as 


where  X..  is  the  location  of  point  A  in  the  ith  dimension  and  X.„  is  the 
lA  —  iB 

location  of  point  B  in  the  ith  dimension. 


This  rationale  can  be  used  to  classify  some  point,  T,  as  being  the 


member  of  one  of  n  classes  (A^  ,  j  *  l,n),  by  finding  the  minimum  vector 


length  from  T  to  A^  (j»l,n).  In  other  words  T  is  said  to  be  a  member  of 


the  class  which  is  closest  to  it,  on  the  average,  across  all  p  dimensions. 
The  minimum  vector  length  is  defined  as 


Notice  that  equation  3  has  no  provision  for  variability  in  the  n  classes. 


therefore,  L  .  is  chosen  as  being  the  shortest  linear  vector  length.  If 
min 

each  class  has  the  sane  variability  associated  with  it  this  causes  no 
difficulty.  In  this  experiment,  however,  there  were  considerable  differ¬ 
ences  in  variability  between  the  classes  so  that  a  modification  of  equation 
3  had  to  be  made.  The  standard  deviation  (SD)  was  used  to  modify  the 
distance  between  T  and  related  to  each  dimension  thus  removing  the 
effects  of  variability  from  each  class.  This  normalized  minimum  distance 
equation  is  expressed  as 


min’ 


min  [  £ 

j-l,...,n  Li-1  L 


xn  -  W/SDij 


j!] 


(4) 


where  SD^  is  the  standard  deviation  associated  with  the  j_th  class  in  the 
ith  dimension. 


In  the  application  of  this  method  to  the  classification  of  mineralogy, 
the  "dimensions"  are  spectral  bands  or  ratios  of  spectral  bands  and  the 
"classes"  are  mineralogical  types.  Eight  spectral  bands  (Table  3)  and  all 
possible  unique  ratios  of  those  spectral  bands  were  used  to  classify  the 
mineralogical  type  of  an  input  sand  as  one  of  five  categories.  The  object 
was  to  make  each  category  as  homogeneous  as  possible  so  that  the  moisture 
and  grain  size  regressions  which  followed  would  be  sensitive  to  small  scale 
spectral  changes. 


The  Prediction  of  Moisture  and  Grain  Size 

Using  the  AQUASAND  model  we  found  that  information  related  to  moisture 
content  of  sands  is  best  derived  from  the  spectral  region  beyond  1.0  urn. 
This  is  the  result  of  the  spectral  reflectance  of  sand  in  this  region  being 
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reduced  by  absorption  in  proportion  to  the  amount  of  water  present. 
(Exceptionally  high  spectral  absorption  is  noted  near  1.4  and  1.9  pm.) 
Although  the  spectral  reflectance  in  these  regions  is  highly  correlated  to 
moisture  we  did  not  consider  them  since  atmospheric  absorption  prohibits 
their  use  by  an  airborne  sensor. 

Changes  in  grain  size  seem  to  manifest  themselves  most  clearly  in  the 
shorter  wavelengths  (0.4-0. 7  pm).  Grain  size  information  is  gained  by 
light  being  reflected  from  sand  grains  below  the  surface  through  surface 
grains.  The  transmittance  through  the  surface  grains  is  reduced  by 
internal  scattering  and  absorption  of  the  particle.  Since  both  of  these 
factors  are  dependent  on  thickness,  the  bulk  reflectance  of  a  sand  is 
dependent  to  some  degree  on  the  grain  size.  Theoretically  a  coarsely  grained 
sand  should  have  a  lower  reflectance  than  a  fine  grained  sand  of  the  same 
mineralogical  composition,  surface  frosting,  and  with  similar  moisture 
content.  According  to  our  measurements  this  appears  to  be  the  case. 

This  grain  size  phenomenon  can  be  confounded  in  two  ways.  First,  if 
there  is  no  scattering  or  absorption  within  the  grains  (i.e,,  a  perfectly 
clear  material  at  all  wavelengths)  there  can  be  no  attenuation.  Fortunately, 
even  in  our  purest  quartz  sands  there  were  enough  impurities  and  inclusions 
to  give  some  attenuation.  Second,  the  sand  grains  may  be  opaque  and  thus 
attenuate  too  much  light.  This  appears  to  be  the  case  in  the  heavy  mineral 
and  carbonate  beaches.  Most  of  the  bulk  reflectance  for  these  two  types 
was  the  result  of  surface  refectance  and  essentially  none  from  light  trans¬ 
mitted  through  the  surface  grains  from  below.  We  were  unable  to  create 
accurate  grain  size  equations  for  these  types. 


Ucilizing  Che  physical  phenomena  discussed  above  we  were  able  Co 
develop  mulciple  linear  regression  equacions  for  prediccing  moisCure  in 
all  five  mineralogical  classes  and  grain  size  for  chree  of  Che  five 
mineralogical  classes.  The  basis  for  all  Che  regressions,  except  one, 
was  Che  sample  group  corresponding  Co  a  given  mineralogical  class.  The 
single  excepcion  was  Che  grain  size  equaCion  corresponding  Co  a  pure  quarcz 
beach.  Our  samples  wichin  Chis  Cype  consisted  of  a  single  grain  size 
(0.22  nnu)  and,  as  such  did  noC  provide  an  adequace  basis  for  regression 
equacions.  For  chis  case  we  used  AQUASAND  generaced  spectra  Co  simulate  a 
wide  range  of  grain  sizes  in  order  to  add  grain  size  variability  to  the 
data  sec. 

Seventeen  spectral  bands  becween  0.4  and  2.5  ym  were  chosen  for  use  in 

the  regressions  (Table  2).  Within  the  17  bands,  only  those  which  were 

predicted  by  the  AQUASAND  model  to  be  most  informative  were  used.  In  this 

way  we  could  be  reasonably  certain  that  the  regression  equations  would 

respond  the  correct  parameter  and  thus  yield  accurate  predictions.  The 

predictive  equations  together  with  the  associated  standard  errors  (SE) ,  and 

2 

coefficient  of  variation  (R  )  are  given  in  Table  4. 

In  summary,  the  MOGS  algorithm  represents  (see  Figure  1)  a  computer 
controlled  package  of  equations.  The  input  is  a  set  of  17  spectral 
reflectance  bands  obtained  from  an  unknown  sand.  Based  on  these  bands, 
the  sand  is  classified  as  being  a  member  one  of  five  mineralogical  types. 
Depending  on  the  mineralogical  type,  the  appropriate  moisture  and  grain 
size  (where  applicable)  equations  are  applied  to  the  data.  The  output  from 
the  MOGS  algorithm  is  the  predicted  mineralogical  class,  the  predicted 
moisture,  and  the  predicted  grain  size. 
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TABLE  4 


Multiple  linear  regression  equations  for  the 
prediction  of  moisture  and  grain  size.  The 
equations  are  listed  by  mineralogical  class. 
Grain  size  is  in  mm. 


a.  Iron  stained  quartz  -  Atlantic  coast 
Predicted  moisture  7.  ■  67.964  -  65.046  (gyjjj-^41) 

S.E.  -  ±  3.082,  R2  -  0.888 

Predicted  grain  size  ■  6.37  -  3.4634  (Band  7)*,/,il  +  .0300  (Band  1)  +  .01672(Band  15) 
2 

S.E.  “  +  0.13  nun,  R  »  0.603 

b.  Iron  stained  quart2  -  Michigan  coast 

Predicted  moisture  7.  -  60.149-49.961  (§“rxf> *2 ■ 226  (Iand  l'^ 

S.E.  -  ±  2.562,  R2  -  0.970 


Predicted  grain  size  -  0.6405-0.0152  (Band  5)  -  .0047  (Band  17) 

2 

S.E.  -  +  0.055  mm,  R  -  0.558 

c.  don- Iron  stained  quartz 

Predicted  Moisture  7.  -  127.02-65.159  (g*"j  -  64.054  (|^  ^) 

S.E.  -  +  2.122,  R2  -  0.971 

Predicted  grain  size  -  1.158-2.328  (Band  10)  +  .3201  (f^j-y)*  0.2C58  (Band  10 
S.E.  and  R^  not  applicable 

d.  Carbonate 

Predicted  moisture  1  -  596.28-642  [■)-  1.081  (Band  14)  +  0.1538  (Band  17) 

S  E.  ■  +  4.092,  R2  r  0.879  i  f>o  grain  size  equation. 


e.  Heavy  mineral 

Predicted  moisture  2  -  19.284  +  11.194  (|^j  yy)  -  1.081  (Band  14)  +  0.1538  (Band  17) 

S.E.  -  +  4.092,  R2  ■  0.879  No  grain  size  equation 

*  Iron  stain  quartz  was  determined  by  visual  examination  under  stereo-microscooe. 
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TEST  RESULTS  OF  THE  MDGS  ON  LABORATORY  SPECTRA 
The  MDGS  algorithm  was  first  tested  on  70  of  the  81  samples  frcm  which  it 
was  derived  and  the  results  were  very  promising.  Eleven  samples  were  discarded 
because  their  grain  sizes  were  abnormally  large  or  percent  moisture  was  inaccurate. 
The  classification  of  mineralogy  was  99+7.  correct.  The  overall  correlation  of 
predicted  to  actual  moisture  was  967.  (significant  at  the  0.001  level)  and  the 
overall  correlation  of  predicted  to  actual  grain  size  was  887.  (significant  at 
the  0.001  level).  However,  testing  any  equation  or  algorithm  cn  the  samples 
frcm  which  it  was  derived  is  not  conclusive.  For  this  reason  the  MDGS  algorithm 
was  tested  on  several  other  beach  sand  samples  which  were  independently 
collected  and  spectrally  measured  following  the  algorithm  construction. 

These  results  are  given  in  Table  5.  In  each  case  the  MOGS  algorithm  selected 
a  mineralogy  which  allowed  the  moisture  and  grain  size  regression  to  operate 
correctly.  The  independent  test  yielded  an  actual  moisture  to  predicted 
moisture  correlation  of  0.95  (significant  at  the  0.01  level).  The  prediction 
of  grain  size  was  in  no  case  more  than  0.07  ram  different  from  the  actual 
grain  size. 

TEST  OF  THE  MOGS  ALGORITHM  ON  MSS  DATA 

The  next  logical  test  for  the  MOGS  algorithm  was  to  evaluate  it  on 
actual  multispectral  scanner  (MSS)  data.  Such  an  investigation  was  carried 
out  using  aircraft  data  obtained  from  the  ERIM-MSS.  The  ERIM-MSS  (Hasell, 
er.  al,  1974)  is  an  optical  scanning  device,  i.e.,  it  receives  signals 
continuously  while  the  sensor  instantaneous  field  of  view  (IFOV)  moves  over 
the  flight  path  scanning  the  terrain.  The  MSS  is  a  passive  device,  meaning 
it  senses  energy  originating  from  the  sun  and  reflected  fronj  the  terrain. 

The  data  from  the  MSS  used  in  this  analysis  covered  12  discrete  spectral 
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TABLE  5 


Comparisons  of  actual  parameters  to  predicted 
classifications  by  the  MOGS  algorithm.  Samples 
used  here  were  collected  independently  of  those 
on  which  the  algorithm  is  based. 


Sample 

Mineralogy 

Moisture  7, 

Mean 

Grain  Size 

(mm) 

Actual 

Predicted 

Actual 

Predicted 

Percent 

Different 

Actual 

Predicted 

mm 

Different 

MICH1 

MICH 

37.7 

25.6 

-12.1 

.26 

.28 

+  .02 

MICH2 

MICH 

3.2 

3.7 

+  .5 

.28 

.30 

+  .02 

MICH3 

MICH 

0.3 

0.0 

-  .3 

.25 

.32 

+  .07 

MICK4 

MICH 

12.1 

16.0 

+  3.9 

.23 

.29 

+  .06 

MICH5 

MICH 

15.0 

12.7 

-  2.3 

.25 

.30 

+  .05 

MICH6 

MICH 

28.0 

26.2 

-  1.8 

.28 

.29 

+  .01 

TABLE  6 

Spectral  bands  used  by  ERIM-MSS  for 
testing  the  MOGS  algorithm 


Band 


mber 

Range  (pm) 

1 

0.40-0.44 

2 

0.43-0.46 

3 

0.45-0.49 

4 

0.48-0.53 

5 

0.51-0.57 

6 

0.54-0.63 

7 

0.60-0.72 

8 

0.66-0.84 

9 

0.78-1.1 

10 

1.2  -1.4 

11 

1.5  -1.8 

12 

'2.0  -2.6 
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bands  in  the  0.4  to  2.5  ym  region  of  the  electromagnetic  spectrum  (see 
Table  6). 

Field  Test 

Multi-Spectral  Scanner  data  were  collected  over  three  test  sites  along 
the  Michigan  shoreline;  Muskegon  State  Park,  Pentwater  State  Park,  and  Mason- 
Oceana  (M-O)  County  line.  The  flight,  part  of  an  ERIM  flight  test  for  a  NASA 
and  NOAA  program,  took  place  on  1  November  1978  under  clear  skies  between 
1:00  and  2:00  P.M. 

Coincident  with  the  fly-over,  the  following  ground  truth  survey  was 
carried  out  at  the  Pentwater  test  site: 

1.  sand  samples  were  collected  to  confirm  grain  size  and 
moisture  conditions  t\t  the  time  of  flight,  and 

2.  reflectance  panels  (black  and  gray)  were  placed  on  the 
beach  and  radiance  and  irradiance  spectral  measurements 
obtained. 

Additionally,  three  sand  samples  were  taken  at  the  M-0  County  line  site 
and  two  sand  samples  obtained  along  with  a  reflectance  measurement  at  the 
Muskegon  State  Park  test  site. 

The  Michigan  beach  test  was  the  first  test  of  the  MOGS  algorithm  on 

actual  MSS  data.  The  bands  available  on  the  ERD1-MSS  (M7)  as  it  presently 

functions  necessitated  a  change  in  the  band  classification  using  the  MOGS 

algorithm.  The  algorithm  was  changed  to  handle  a  total  of  12  instead  of  17 
bands.  This  was  simply  done  by  truncating  the  five  bands  from  the  algorithm 
not  available  on  the  ERIM-MSS(M7) .  The  resulting  predictions  of  mineralogy 
using  12  instead  of  17  were  not  affected.  Us'ry  only  five  bands  and  ten  ratios, 
the  algorithm  classified  the  mineralogies  very  well.  Tbi>  indicates  that  the 
use  of  eight  bands  and  26  ratios,  as  in 
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Che  original  algorithm,  is  probably  over  precision.  In  addition,  the  117 
bands  give  very  nearly  the  same  predictions  for  moisture  and  grain  size  for 
the  Michigan  beaches  as  the  17  bands  which  were  used  to  develop  the  original 
equations. 

Aircraft  MSS  Data  Reduction 

Table  6  lists  the  12  channels  of  MSS  data  used  in  the  Michigan  tests. 
These  channels  (as  mentioned  previously)  were  selected  as  being  very  close 
to  the  channels  used  in  the  development  of  the  MOGS  algorithm.  The  nine 
channels  corresponding  to  the  shorter  wavelengths  utilized  a  photo-multiplier 
detection  system  while  the  three  longer  wavelength  channels  used  an  InSb 
(lead  sulfide)  detector.  The  data  were  recorded  digitally  on  a  High  Density 
Digital  Tape  (HDDT)  for  later  computer  processing. 

Using  an  ERIM  computer  video  display,  the  exact  location  of  the  test 
sites  were  determined  and  the  data  for  all  three  test  sites  were  transferred 
to  a  computer  compatible  tape  (CCT)  for  use  on  the  University  of  Michigan 
AMDAHL  470  computer  system  (MTS).  Rather  than  converting  all  the  data,  only 
an  area  100  pixels  on  either  side  of  the  test  site  was  recorded  for  each 
test  site.  Each  pixel  for  the  test  data  (600  meter  altitude)  is  approxi¬ 
mately  a  1.5  meter  square.  The  use  of  only  100  pixels  on  either  side  of  the 
test  site  negated  having  to  correct  the  aircraft  data  for  scanning  angle 
effects  of  the  sensor. 

Using  the  generated  CCT  on  MTS,  gray  maps  (computer  generated  images) 
•..ere  produced  to  ascertain  t'..e  exact  areas  from  which  the  sand  samples  were 
collected.  Upon  specific  location  of  the  sample  areas  within  the  test  site, 
the  actual  digital  values  obtained  by  the  remote  sensor  were  extracted.  The 
data  values  obtained  from  MTS  corresponding  to  the  sample  areas  were  then 
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tabulated  and  "calibrated".  This  was  accomplished  by  using  an  MTS  module 
to  extract  the  calibration  values  related  to  each  scan  line  which  contained 
a  sample  point.  The  "black  body"  calibration  value  was  subtracted  from  the 
sample  value  for  each  sample  to  normalize  the  data.  We  found  essentially 
no  variability  between  scan  lines  in  terms  of  the  calibration  values. 

Concurrent  with  the  data  extraction,  the  sand  samples  obtained  from 
the  test  site  were  measured  for  reflectance  using  the  Cary  14  spectrophoto¬ 
meter.  The  signal  values  obtained  from  the  scanner  data  were  then  plotted 
against  the  appropriate  band  reflectance,  as  measured  with  the  Cary  14,  to 
achieve  a  signal  value  to  reflectance  transfer  curve.  For  some  reason, 
probably  related  to  the  non-Lamb ertian  nature  of  objects  in  the  scene,  the 
relationship  did  not  appear  to  be  strictly  linear  in  all  bands.  This  non- 
lipearity  appears  to  be  most  pronounced  for  the  near  infrared  spectral 
bands  which  use  the  InSb  detector.  In  addition,  the  reflectance  calibration 
panels  which  were  deployed  at  the  Pentwater  Site  do  not  plot  well  with  the 
sand  values.  This  means  that,  in  this  case,  the  panels  were  not  particularly 
useful  in  calibrating  the  MSS  system  to  values  of  reflectance.  Our  final 
signal-value-to-surface  reflectance  calibration  was  accomplished  using  the 
reflectance  of  a  few  sand  samples  (four)  plotted  against  the  MSS  digital 
signal  value.  The  digital  values  of  other  areas  on  the  beach  were  then 

altered  to  reflectance  using  these  curves.  Thus  the  ERIM-MSS (M7)  received 
radiance  value  was  converted  to  a  reflectance  value  using  this  technique. 

Aircraft  Test  Results 

The  test  results  using  the  aircraft  sinner  data  are  encouraging 
(Table  7).  The  correlation  (R)  of  predicted  to  actual  moisture  content 
is  0.91  (significant  at  the  0.01  level)  and  the  prediction  of  grain  size  is, 
in  no  case,  greater  than  0.09  mm  different  from  the  actual  grain  size.  The 


18 


TABLE  7 


Comparison  of  acCual  parameters  to 
predicted  classifications  by  the 
MOGS  algorithm.  The  sand  spectra 
used  here  were  collected  by  the  ERIM-MSS 


Sample 

Mineralogy 

Moisture  X 

Grain  Size  (mm) 

Percent 

mm 

Actual 

Predicted 

Actual 

Predicted 

Different 

Actual 

Predicted 

Different 

Pentwater 

MICH 

22.3 

12.1 

-10.2 

.25 

.32 

+ 

.07 

Pentwater 

MICH 

1.0 

o:o 

-  1.0 

.23 

.28 

+ 

.05 

Pentwater 

MICH 

8.0 

10.1 

+  1.9 

.22 

.25 

+ 

.03 

Pentwater 

MICH 

28.0 

16.9 

-11.1 

.26 

.35 

+ 

.09 

Pentwater 

MICH 

5.0 

2.2 

-  2.3 

.24 

.26 

+ 

.02 

M-0  Line 

MICH 

Not  Measured 

.29 

.32 

+ 

.06 

M-0  Line 

MICH 

Not  Measured 

.27 

.36 

+ 

.09 

>1-0  Line 

MICH 

Not  Measured 

.29 

.37 

+ 

.08 

Muskegon  St. 

Pk.  MICH 

Not  Measured 

.25 

.30 

+ 

.05 

predicted  grain  sizes  shown  in  Table  7  are  all  larger  than  the  actual 
measured  grain  size.  The  statistical  analysis  on  the  sand  samples  used  to 
create  the  Michigan  grain  size  regression  equation  showed  the  sands  to  be 

moderate  to  well  sorted  (i.e. ,  low  standard  deviation  of  grain  size)  but 
exhibiting  a  negative  skewness  in  almost  all  cases.  This  negative  skewness, 
indicative  of  a  coarse  fraction,  could  have  biased  the  grain  size  prediction 
algorithm  and  caused  the  larger  than  actual  prediction  of  grain  size  to 
result. 

The  moisture  prediction  appears  particularly  poor  for  high  moisture 
contents.  This  is  most  likely  the  result  of  wet  sand  at  the  test  sites 
exhibiting  bi-directional  dependencies  (i.e.,  a  failure  to  behave  in  a 
Lambertian  manner).  •  These  bi-directional  characteristics  are  enhanced  at 
low  sun  angles  and  are  not  accounted  for  by  the  MOGS  algorithm.  Although 
the  flight  Look  place  at  1:30  EST  the  sun  was  only  38°  above  the  horizon  on 
November  1. 

A  laboratory  experiment  was  conducted  to  quantify  the  effects  of  bi¬ 
directional  dependencies  of  wet  sand.  Figure  2  shows  the  laboratory  set-up 
used  to  conduct  the  test.  The  sand  was  spread  evenly  in  a  70  x  70  cm  tray 
to  a  depth  of  2.5  cm.  Illumination  was  provided  by  a  200  yw  projector  lamp 
whose  beam  was  optically  collimated  to  provide  a  highly  directional  source(i.e.  =G9) 
The  detector  used  was  a  Coherent  Optics  power  meter  which  utilized  a  silicon 
detector.  The  aperture  of  the  detector  was  reduced  using  baffles  such  that 
an  area  approximately  5  cm  square  on  the  sand  surface  was  sensed.  To  avoid 
having  the  detector  head  shadowing  the  sand  it  was  offset  such  that  it  was 
viewing  at  an  angle  of  23°  from  the  vertical.  As  is  shown  in  Figure  2,  the 
lump  was  moved  through  several  angles  from  normal  to  the  sand  surface  to 

ft 

nearly  parallel  with  it.  At  each  angle  the  flux  to  the  detector  was  measured 

and  recorded.  - 
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According  to  theory,  if  a  surface  is  perfectly  Lambertian  then  the 
exitance  from  the  surface  should  decrease  with  the  cosine  of  the  angle  of 
incidence.  Therefore  the  exitance  observed  at  any  illumination  angle 
should  be  equal  to  the  exitance  with  the  illumination  normal  to  the  surface 
multiplied  by  the  cosine  of  the  angle  of  interest.  Using  uhis  relationship, 
a  useful  mer lure  of  the  Lambertian  nature  of  a  material  is 


R 


Mq  cos  9 


(5) 


where  R  is  the  measure  of  Lambertian  nature,  M0  is  the  exitance  observed 
at  illumination  angle  9  and  MQ  is  the  exitance  observed  with  the  illumi¬ 
nation  normal  to  the  surface.  For  a  perfectly  Lambertian  material  R 
should  be  1.0  for  all  9.  If  R  is  less  than  1.0  then  less  flux  is  being 
reflected  at  angle  9  than  was  expected,  and  if  R  is  greater  than  1.0  the 
reverse  is  true. 

Using  the  experimental  design  just  described,  two  sand  conditions  were 
investigated:  dry  sand  and  sand  with '■>-25%  moisture.  As  can  be  seen  in 

Figure  3A,  dry  sand  is  quite  Lambertian  in  character  at  angles  up  to  60° 
incident  (30°  above  the  horizon).  The  sun  angle  during  the  Michigan  test 
flight  was  approximately  52°  incident  indicating  the  exitance  from  dry  sand 
was  ~84%  of  that  expected  from  a  Lambertian  material.  This  is  in  fact  quite 
good  for  this  incidence  angle  although  a  more  nearly  normal  incidence  angle 
(higher  sun)  would  have  been  desirable. 

The  wet  sand  (Figure  3B)  exhibits  strong  bi-directional  properties 
indicative  of  a  semi-specular  surface.  The  raw  data  indicates  a  very  high 
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FIGURE  3B.  WET  SAND  CASE 


exitance  at  0°  incidence  followed  by  a  rapid  fall  off  ac  other  incidence 
angles  relative  to  that  observed  for  dry  sand.  The  implication  is  that 
there  is  a  specular  return  at  0°  incident — a  very  non-Lambertian  charac¬ 
teristic;  for  such  a  case  the  assumption  of  a  Lambertian  surface  does  not 
hold.  Note  from  Figure  3  that  for  sun  angle  +15®  from  zenith  that  wet 

and  dry  sands  decrease  in  reflectance  approximately  the  same  amount  and 
therefore  minimize  non-Lambertian  effects.  It  should  be  noted  that  the 

goniometric  study  described  was  not  a  spectral  study  but  rather  used  a  light 

source  that -ranges  the  visible  and  is  weighted  to  the  infrared. 

In  summary,  the  Cary  14  Spectrophotometer  makes  a  reflectance  measure¬ 
ment  that  is  hemispherical,  i.e.,  the  reflectance  properties  whether 
bi-directional,  Lambertian,  or  diffuse  have  little  bearing  on  the  reflectance 
that  is  sensed.  Conversely,  the  MSS  only  intercepts  radiation  which  is 
reflected  in  a  certain  direction  as  opposed  to  the  diffuse  sensing  of  a 
spectrophotometer.  For  an  MSS,  this  sensing  is  normal  to  the  sensed  surface 
at  the  center  of  the  scan  and  progressively  more  oblique  toward  the  edges  of 
the  scan.  The  type  of  illumination  is  highly  variable  ranging  from  diffuse 
with  an  overcast  sky  to  nearly  specular  on  a  clear  sunny  'day.  In  the  former 
case,  the  geometry  is  much  like  that  of  a  spectrophotometer  in-  reverse — 
diffuse  source  and  specular  sensing.  Because  of  reciprocity  the  apparent 
reflectance  as  derived  from  MSS  data  ignoring  path  radiance  factors  should 
be  exactly  the  same  as  the  reflectance  sensed  by  a  spectrophotometer.  Thus, 
to  minimize  bi-directional  reflectance,  future  aircraft  flights  should  be 
made  during  complete  mid-altitude  (3000  m)  cloud  cover  or  sunny  skies  with 
the  sun  close  to  the  zenith  (summer  sun)  '.a  indicated,  from  Figure  3. 

Applying  the  same  moisture  and  grain  size  equations  used  in  the  previous  ■ 
analysis,  the  entire  Pcntwater  State  Park  beach  on  a  pixal-by~pixal  baait  (in 
this  case  1,5  x  1.5  meters)  was  classified  with  respect  to  grain  size  and 
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moisture  content:.  The  two  MOGS  generated  digital  maps  (see  Figure  4)  show 
the  predicted  moisture  and  grain  size  distribution  on  the  beach  at  Pentwater. 

Also  included  on  Figure  4  is  a  panchromatic  aerial  photograph  take  coincident 
with  the  MSS  data.  The  ground  truth  measurements  taken  at  the  time  of  flight 
correlate  well  with  these  images.  Figure  4  helps  to  demonstrate  how  potentially  an 
entire  sandy  coastline  could  be  analyzed  in  respect  to  moisture,  grain  size, 
and  gross  mineralogy  using  a  small  subsection  as  calibration. 

CONCLUSION 

The  development  of  the  MOGS  algorithm  has  demonstrated  the  feasibility 
of  obtaining  quantitative  moisture  and  grain  size  information  from  the 
spectral  reflectance  of  beach  sands.  The  determination  of  grain  size  is 
dependent  on  the  sand  grains  being  neither  opaque  nor  perfectly  clear. 

The  two  stage  nature  of  the  MOGS  algorithm  is  directly  responsible  for 
its  broad  applicability  without  loss  of  detail.  By  separating  the  minera- 
logical  types  prior  to  the  prediction  of  moisture  and  grain  size  much  of 
the  variability  which  could  easily  hide  small  scale  changes  is  removed. 

The  use  of  a  vector-length  discriminant  function  to  classify  mineralogy 

worked  extremely  well  in  this  application,  since  36  different  dimensions 
could  be  evaluated  simultaneously.  A  multistage  approach  involving  multiple 

classification  techniques  is  a  powerful  tool,  one  which  may  be  very  useful 
in  many  areas  of  remote  sensing.  _ _ 

The  Lake  Michigan  field  test  has  further  demonstrated  the  MOGS  algorithm’s 
applicability  to  remotely  sensed  field  data.  Grain  size  was  predicted  to 
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MOISTURE  DISTRIBUTION  IMAGE  -  LIGHTER  SHADE  INDICATES  HIGHER  MOISTURE 
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GRAIN  SIZE  DISTRIBUTION  IMAGE  -  LIGHTER  SHADE  INDICATES  LARGER  GRAIN  SIZE 


PANCHROMATIC  AERIAL  PHOTOGRAPH 


5 JOITAL  IMAGERY,  GENERATED  BY  THE  MOGS  ALGORITHM  SHOWING  THE  DISTRIBUTION 
OF  MOISTURE  AND  GRAIN  SIZE  ON  PENTWATER  BEACH  (PENTWATER  STATE  PARK 
MICHIGAN)  .  WHITE  AREAS  ARE  EITHER  UNCLASSIFIED  REGIOHS  OR  OPEN  WATER. 


within  0.09  mra  mean  diameter  of  the  actual  size  while  beach  moistures  below 
20%  were. accurately  predicted.  In  all  cases  the  computer  algorithm  correctly 
identified  the  Michigan  beach  mineralogy  as  being  a  predominately  iron- 
stained,  quartz-feldspar  beach. 
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